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Abstract 

This  thesis  addresses  the  application  issues  raised  implementing  flight  control  designs  derived  from 
optimal  control  theory  and  the  challenges  in  obtaining  acceptable  handling  qualities  when  using  these 
techniques.  Using  the  USAF  TPS  PCS  project  as  the  controller  architecture,  four  controllers  were  designed 
using  classical  methods,  and  H2,  ft-,  and  mixed  H2/ft.  optimal  control  theory.  These  designs  were 
implemented  in  the  Calspan  VSS  11  Learjet,  simulating  unstable  aircraft  longitudinal  dynamics  and  a  limited 
handling  qualities  flight  test  evaluation  was  performed.  The  design  phase  found  the  optimal  control 
techniques,  as  applied,  difficult  to  design  to  handling  qualities  specifications.  The  H2,  H.-,  and  mixed  H2/Ho= 
controllers  were  unstable  and  often  contained  high  frequency  poles,  which  were  difficult  to  implement.  The 
flight  test  rated  the  designs  acceptable  on  approach,  but  no  handling  qualities  level  for  approach  was 
determined. 
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MODERN  FLIGHT  CONTROL  DESIGN,  IMPLEMENTATION  AND  FLIGHT  TEST 


/.  Introduction 

1.1.  Background/Motivation 

Flight  controller  designs  for  piloted  aircraft  have  become  increasingly  more  complex  in  the  past  few 
years  with  the  advent  of  fly-by-wire  control  systems  [Ber88].  Modern  aircraft  now  have  the  capability  to 
use  both  conventional  and  innovative  control  surfaces  and  systems  to  increase  performance  and  reduce  drag 
and  radar  cross  section.  While  these  systems  greatly  enhance  capabilities,  the  challenges  to  the  flight 
control  designer  have  grown  exponentially.  Single-input-single-output  (SISO)  systems  become  multiple- 
input-multiple-output  (MIMO)  systems.  Modeling  the  aircraft  dynamics,  which  has  never  been  an  exact 
science  [DFT92],  is  more  complex.  To  simplify  the  design  process,  a  linear,  time  invariant  system  model  is 
developed  about  an  operating  point  or  nominal  condition.  The  primary  purpose  of  a  flight  control  system 
for  piloted  aircraft  is  to  provide  acceptable  handling  qualities  to  accomplish  a  mission  flight  phase  [TPS95]. 
Embedded  in  this  requirement  for  good  handling  qualities  is  the  objective  to  obtain  nominal  and  robust 
stability  and  performance.  Ultimately  the  pilot  completes  the  flight  control  system,  and  as  handling 
qualities  degrade,  the  pilot  becomes  unable  to  perform  piloting  duties.  Stability  and  performance  alone  are 
not  sufficient  if  a  pilot  cannot  control  the  aircraft. 

There  has  been  considerable  research  developing  design  techniques  to  optimize  flight  control 
systems  to  reduce  noise  introduced  into  the  system,  and  provide  robust  stability  to  the  system  [Gan86]. 
These  techniques  include  H2  control  theory,  H,.  control  theory,  and  a  Mixed  H2/EL  control  theory  [RW95]. 
H2  optimization  minimizes  the  system’s  output  energy  when  the  system  is  faced  with  white  Gaussian  noise 
inputs.  The  H2  design  technique  targets  noise  rejection  as  its  main  objective.  PL,  optimization  minimizes 
the  system’s  output  energy  to  unknown,  bounded  energy  inputs,  which  results  in  a  highly  robust  system. 

The  design  technique  targets  system  stability  margins  and  also  provides  good  tracking.  Mixed  H2/PP~ 
optimization  combines  the  objectives  of  H2  and  PL  optimization  by  limiting  the  II2  problem  with  an  PL 
constraint.  The  Mixed  PI2/PL  design  technique  optimally  trades  off  noise  rejection  with  the  tracking  and 
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robust  stability  objectives  to  obtain  a  system  targeting  all  of  these  goals.  All  of  these  design  methods  are 
capable  of  handling  both  SISO  and  MIMO  systems.  These  techniques  do  not  directly  provide  assurances  of 
acceptable  handling  qualities  in  the  flight  control  systems  designed  by  these  methods.  The  previous 
research  on  these  methods  has  focused  on  improving  the  numeric  solutions  of  these  techniques.  Although 
example  applications  of  these  techniques  have  been  analyzed  for  robustness  and  some  measures  of 
performance,  little  has  been  done  to  investigate  the  resulting  design’s  handling  qualities. 

There  are  several  handling  qualities  design  criteria  for  evaluating  flight  control  designs  described  in 
the  Department  of  Defense  MIL-STD-1797A,  Flying  Qualities  of  Piloted  Aircraft  [Mil90].  These  criteria 
are  used  to  predict  an  aircraft’s  handling  qualities  based  on  historical  pilot  evaluations.  Of  these  criteria, 
Hoh’s  Bandwidth  Criteria  [HohSl]  and  Ralph  Smith  Criteria  [SG79]  are  accepted  as  predictors  of  an 
aircraft’s  handling  qualities.  The  USAF  Test  Pilot  School  (TPS),  as  part  of  its  curriculum,  uses  these 
criteria  to  evaluate  student  designed  flight  control  systems  during  the  Flight  Control  System  (FCS)  design 
project.  This  project  then  implements  the  FCS  in  the  CALSPAN  Variable  Stability  Simulator  11  (VSS  II) 
Learjet  [Bal93]  to  allow  the  TPS  students  to  evaluate  the  handling  qualities  of  their  design.  This  project  is 
well  suited  for  implementation  and  a  limited  handling  qualities  evaluation  of  augmented  longitudinal  flight 
control  systems  designed  with  optimization  techniques. 

1.2  Research  Objectives 

The  purpose  of  this  thesis  is  to  investigate  the  implementation  of  flight  control  systems  designed 
using  H2,  IL,,  and  Mixed  H2/Ii-  design  techniques  and  evaluate  the  resulting  handling  qualities  of  these 
designs.  This  research  applies  modern  optimal  control  theory  to  the  USAF  TPS  FCS  design  project,  used  in 
the  TPS  curriculum  from  1994  through  1995.  These  designs,  as  well  ^  a  classically  designed  flight  control 
system,  were  implemented  in  the  CALSPAN  VSS  n  Learjet  and  evaluated  for  handling  qualities  during 
flight  test.  The  classically  designed  flight  control  system  served  as  a  control  case  for  comparison  and  was 
used  to  help  verify  proper  implementation  of  the  un-augmented  aerodynamics  simulated  by  the  VSS  II 
system.  The  classical  design  serves  as  a  base-line  and  is  not  intended  to  be  compared  to  the  optimal  control 
theory  designs  to  determine  which  method  can  achieve  a  “better”  controller.  This  thesis  addresses  the 
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application  issues  raised  implementing  flight  control  designs  derived  from  optimal  control  theory  and  the 
challenges  in  obtaining  acceptable  handling  qualities  when  using  these  techniques. 

1.3  Thesis  Outline 

This  thesis  is  divided  into  six  chapters,  including  this  introductory  chapter.  The  second  chapter 
presents  the  necessary  groundwork  for  this  research.  The  control  theory  used  to  develop  the  flight  control 
designs  is  discussed  in  the  first  section  of  Chapter  H.  An  explanation  of  the  handling  qualities  criteria  used 
to  evaluate  the  designs  is  then  presented.  The  third  section  of  this  chapter  addresses  the  TPS  PCS  project 
on  which  this  research  is  based.  Finally,  a  description  of  the  CALSPAN  VSS  II  Learjet  and  its  limitations 
is  presented. 

In  Chapter  III,  the  various  models  used  in  this  research  are  developed.  The  aircraft  aerodynamic 
model,  in  both  transfer  function  and  state  space  form,  is  derived  from  the  stability  derivatives  of  the 
intended  target  aircraft.  The  models  for  wind  gust  noise  and  sensor  noise,  used  in  the  design  and  computer 
evaluation,  are  then  developed.  The  third  section  addresses  the  structure  of  the  closed  loop  system  used  in 
each  of  the  designs.  It  also  explains  the  model  used  in  the  computer  simulations  to  evaluate  the  designs  and 
in  the  handling  qualities  predictions.  The  fourth  section  describes  the  set-up  for  each  design  method, 
comments  on  the  design  process  and  presents  the  resulting  designs. 

Chapter  IV  describes  the  implementation  of  the  flight  control  designs.  It  discusses  the  controller 
model  used  in  implementation  and  the  issues  raised  developing  this  model.  This  chapter  also  presents  the 
results  of  the  ground  and  flight  verification  and  validation  testing. 

The  fifth  chapter  addresses  the  flight  test.  It  describes  the  flight  test  objectives,  and  the  methods  and 
procedures  used  during  the  flight  test.  The  flight  test  results  are  presented  and  then  analyzed.  The  flight 
test  analysis  focuses  on  the  handling  qualities  ratings  determined  during  the  test,  and  on  the  conclusions  that 
can  be  drawn  based  on  these  results. 

The  final  chapter  summarizes  the  findings  of  this  study,  presents  some  conclusions  and  provides 
recommendations  for  areas  of  future  research. 
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II.  Preliminaries 


2.1.  Control  Theory 

The  control  theory  and  the  optimization  techniques  presented  in  this  section  are  all  well  developed 
processes.  This  section  will  summarize  the  basic  principles  and  assumptions  made  by  the  design  methods 
used  for  each  of  the  flight  controllers  examined  in  this  thesis.  A  very  general  description  of  the  classical 
control  theory  used  to  design  the  classical  controller  will  first  be  discussed.  The  remaining  subsections 
present  the  basic  development  of  H2,  H-,  and  Mixed  H2/  IL  optimization. 

2.1.1.  Classical  Control  Theory 

Flight  control  designs  have  classically  used  single  loop  frequency  response  and  root  locus  design 
techniques  to  solve  the  control  problem.  When  multiple  feedback  loops  are  required  to  achieve  the  desired 
stability  and  performance  objectives,  classical  control  theory  predominantly  relies  on  one-loop-at-time 
frequency  response  and  root  locus  design  techniques  [Gan86].  While  classical  control  theory  encompasses 
a  broader  range  of  techniques  , then  presented  here,  this  discussion  will  be  limited  to  the  techniques  specified 
in  the  TPS  project  outlined  later  in  this  chapter. 

The  classical  design  method  starts  with  a  block  diagram  like  that  shown  in  Figure  2.1.  Here  G 
represents  the  linearized  aircraft  modeled  dynamics  and  K  the  designed  controller.  The  commanded  input  is 
labeled  r,  and  the  outputs  of  interest  are  depicted  by  y.  Transfer  function  representations  for  G  and  K  are 
commonly  used,  as  is  a  negative  feedback  sign  convention.  The  controller  can  be  a  simple  gain. 


G 

3L 

ih 

K 

— 

Figure  2.1.  Block  Diagram 


2-1 


or  a  dynamic  filter.  A  root  locus  of  all  possible  closed  loop  root  locations  for  a  given  loop  transfer  function 
as  loop  gain  is  varied  is  then  plotted  for  the  system  and  gains  are  chosen  to  place  the  closed  loop  poles  in  a 
desired  region  of  the  plot.  Frequency  response  plots,  or  Bode  plots  are  used  to  analyze  the  loop  transfer 
function  gain  magnitude  and  phase  angle  to  determine  the  system  loop  shape.  The  feedback  gains  and 
controller  dynamics  are  varied  until  the  desired  system  response  is  obtained.  As  multiple  loop  closures  and 
higher  order  controller  dynamics  are  used,  this  iterative  process  requires  significant  engineering  insight  and 
experience  [Gan86]. 

2.1.2.  H2  Control  Theory 

The  objective  of  an  H2  design  is  to  synthesize  a  controller  which  minimizes  the  energy  of  a  system 
output  faced  with  a  zero-mean.  Unit  intensity  White  Gaussian  Noise  (UWGN)  input.  H2  optimization  is  a 
generalization  of  the  standard  Linear  Quadratic  Gaussian  (LQG)  problem  [Doy89].  The  feedback  system 
used  for  H2  optimization  can  be  depicted  in  general  form,  shown  in  Figure  2.2.  The  linear,  time-invariant 
aircraft  model  is  defined  by  P.  The  block  represented  by  P  also  includes  parameter  weightings  used  to 
obtain  the  desired  H2  controller  design.  The  controller  developed  by  the  design  techniques  is  represented 
by  K2.  The  exogenous  inputs  to  the  system  are  shown  as  w,  and  the  controlled  system  outputs  are 
represented  as  z.  The  measured  outputs  from  the  system  are  shown  as  y  and  are  input  to  the  controller. 

The  controller  output  signal  is  labeled  u  and  is  fed  back  into  the  system.  The  symbol  represents  the 
closed  loop  system. 
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zw 


Figure  2.2.  General  H2  Feedback  System 


The  plant  P  can  be  partitioned  into  separate  input  and  output  transfer  functions 


where 


z  =  P^w+P^^u 
y  =  PywW+PyuU 


A  state  space  realization  of  P  can  be  written  as 

X2  =  A2X2  +  B^W+B^U 
Z  =  C^X2+D^W+D^^U 

y  =  Cy^X2+Dy^W+Dy„U 


(2.1) 


(2.2) 


(2.3) 


where  the  2  subscript  indicates  the  problem  setup  for  H2  optimization.  Minimizing  the  energy  of  the  system 
output  to  a  UWGN  input  is  equivalent  to  minimizing  the  2-norm  of  the  closed-loop  transfer  function, 
[RW95].  This  minimal  value  of  the  2-norm,  a,  is  defined  to  be  [Wal94]: 


a=  inf 

K(,s)Stabilizing 

=  inf 

K(s)StabiUzing 

=  inf 

K(s)Stabilizing 


‘ZW 


P^  +  P^K(l-P^IC)-^P. 


yu- 


yw 


(2.4) 

(2.5) 

(2.6) 
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The  following  assumptions  are  made: 


i.  jDjh,  =  0 

ii.  Dyu  =  0 

iii.  {A,  is  stabilizable  and  (Cj,  A)  is  detectable 

iv.  =  I  and 


A-i(OI  S„  ■ 


has  full  column  rank 


vi. 


A  -  j(0l 

^yw 


has  full  row  rank 


Assumption  (i)  assures  that  the  closed  loop  two-norm  of  the  system  is  finite.  Condition  (ii)  makes  the 
development  easier  but  can  be  removed  completely.  Condition  (iii)  is  necessary  for  the  existence  of 
stabilizing  solutions.  Assumption  (iv)  is  a  regularity  condition  which  insures  that  there  is  a  direct  penalty  on 
all  controls  and  no  perfect  measurements.  Finally,  requirements  (v)  and  (vi)  are  required  to  ensure  the 
existence  of  stabilizing  solutions  to  the  two  Algebraic  Ricatti  Equations  (ARE)  in  the  H2  solution. 

All  stabilizing  H2  controllers  can  be  parametrized  by  a  family  of  lower  fractional  transformations 
(LFT)  of  a  transfer  function  J  and  a  constrained  freedom  parameter  2  e  H2  [Wal94],  shown  in  Figure  2.3. 


Figure  2.3.  H2  System  with  Parametrized  Controller 
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and 

+  -Y^cIc^Y^  +  bJI=Q  (2.16) 

where 

S,  =  fi,(/-D^Z)J  ■■  (2.17) 

The  controller  is  unique  and  optimal  if  (2  =  0,  otherwise  the  set  of  all  controllers  such  that  [|rj^||  <  a  is 
parametrized  by  any  Qe  B.2  such  that 

\\Q\\l<a^-a^  (2.18) 

When  Q  is  chosen  to  be  identically  equal  to  zero,  the  resulting  optimal  controller  is 
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(2.19) 


Kf' 

0 

2.1.3.  Control  Theory 

The  objective  of  an  H„  design  is  to  synthesize  a  stabilizing  controller  which  minimizes  the  maximum 
energy  of  the  output  e,  given  a  bounded  energy  input,  d  [Dec94].  The  feedback  system  used  for  KL 
optimization  is  depicted  in  Figure  2.4  and  is  identical  to  the  H2  block  diagram  with  the  input  and  output  re¬ 
labeled  for  distinction  between  the  design  methods. 


The  plant  P  can  be  partitioned  into  separate  input  and  output  transfer  functions 


(2.20) 


where 


e  =  P,dW+P,^u 
y  =  PydW+  P^u 


A  state  space  realization  of  P  can  be  written  as 

=  A„x„+Bdd  +  B^u 
e  =  +  D,dd  + 

y  =  CyX„+Dydd  +  Dy„u 


(2.21) 


(2.22) 
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where  the  °o  subscript  indicates  the  problem  setup  for  H  „  optimization.  Minimizing  the  energy  of  the 
output  to  an  exogenous,  bounded  energy  input  is  equivalent  to  minimizing  the  -  norm  of  the  closed-loop 
transfer  function,  [RW95].  This  minimum  achievable  value  of  the  <=0  -  norm,  y,  is  defined  to  be 

7  =  inf  sup  H2  =  inf  ll^d  ||  (2.23) 

-*  K Stabilizing  K Stabilizing  “ 

where  the  infinity-norm  of  is 

lkj|L  =  sup  oIt,,]  (2.24) 


and  <7  denotes  the  maximum  singular  value. 
The  following  assumptions  are  made: 

i.  Dgj=  0 

ii.  Dyu  =  0 


iii.  {A,  Bu)  is  stabilizable  and  (Cy,  A)  is  detectable 


iv-  =  I  and  DyjDj.^=I 


V. 


VI. 


A  -  jcol 

A  -  jojl 


D 


yd 


has  full  column  rank 


has  full  row  rank 


Assumption  (i)  and  (ii)  are  not  required  for  a  solution  to  exist,  but  reduce  the  complexity  of  the  solution. 
Condition  (iii)  is  necessary  for  the  existence  of  stabilizing  solutions.  Assumption  (iv)  is  a  regularity 
condition  which  insures  that  there  is  a  direct  penalty  on  all  controls  and  that  no  perfect  measurements  are 
allowed.  Finally,  requirements  (v)  and  (vi)  are  required  to  ensure  the  existence  of  stabilizing  solutions  to 
the  two  AREs  in  the  H,„  solution. 

The  E>.  optimal  controller  found  by  the  optimization  process  is  not  unique.  This  process  is  iterative 
and  is  based  on  the  parametrization  of  all  sub-optimal  controllers  where 

k<t|L<r  (2.25) 

The  family  of  all  admissible  controllers  which  satisfy  (2.25)  is  given  by  the  LFT  (Figure  2.5) 
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Figure  2.5.  IL.  System  with  Parametrized  Controller 


where 


and 


Jis)  = 


uy 


vy 


Jur' 

Kf 

Jyr_ 

= 

-K, 

0 

/ 

I 

0 

Aj  =A-KjC^-  +y-^Y^Cj  (c,  - 


K,  =  [bI 

Kf=Y^Cl+B,Dl^ 

K,,  =  -(r-^D^.,B^X^  +  C,_  ](l-y-h„xX 

Kfl=y-^Y„C^D,,  +  B,^ 

The  matrices  and  K»  are  solutions  to  the  AREs 


(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 


(A-B,^DlC,fX„  +  X^{A-B,_Dj^C,)  +  X„(y-^BjB^-B,jl)x^  +  CjC,  =  0  (2.32) 

where 

Ce={l-D,„Dl)c,  (2.33) 

and 
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where 


(A  -  )Y„  +  Y^(A~  B^Dl^Cy_  f  +  Y„(y-^Cj C,  -  C[_  C,_  )y^  +  bJ^  =  0  (2.34) 

B,  =  B^(l  -  (2.35) 

With  that  stated,  Q  can  be  chosen  to  be  any  2  e  H«.  such  that 

felL  <  r  (2.36) 

This  parametrization  of  a  controller  K  is  valid  if  and  only  if  ±e  following  three  conditions  hold  : 

i.  Hx  e  dom(Ric)  with  =  Ric(Kx)  ^  0 

ii.  Hy  e  dom(Ric)  with  Y^  =  RicQiY)  ^  0 

iii.  p(Y^XJ<y^ 

where  Hx  and  Hy  are  the  Hamiltonians  associated  with  (2.32)  and  (2.34)  [Wal94]. 

As  previously  stated,  the  process  of  finding  a  controller  with  a  closed-loop  infinity-norm  close  to  y  is 
iterative.  A  value  for  j  is  selected,  then  the  above  three  conditions  are  checked.  If  any  one  of  the 
conditions  fail,  increase  y  and  repeat  the  check.  If  all  three  conditions  are  met,  decrease  y  until  the  desired 
accuracy  is  obtained. 

2.1.4.  Mixed  H2/H„  Control  Theory 

The  objective  of  mixed  H2/Ho„  optimization  is  to  achieve  the  robustness  of  H,„  control  with  the  noise 
rejection  properties  of  H2  control.  The  mixed  diagram  is  shown  in  Figure  2.6.  This  method, 
described  by  Ridgely  [Rid91]  and  Walker  [Wal94],  poses  the  H2  problem,  to  minimize  a  with  an  IL 
constraint  that  is  constraining  the  infinity-norm.  The  final  mixed  design  goal  is  to  determine  an  admissible 
controller  that  meets  the  following: 

inf  |r^  I  ,  subject  to  the  constraint  I  <7  (2.37) 
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Figure  2.6.  General  Mixed  H2/H,.  Design  Diagram 


The  mixed  system  can  be  written  in  general  form  as 


e 

'Ped 

P 

*  ew 

Peu 

z 

= 

Pzd 

P 

zw 

Pzu 

w 

y_ 

p 

yw 

- 1 

3 

u 

(2.38) 


A  state  space  representation  is  given  by 

x=  Ax  +  B^d  +  B^^w+B^u  (2.39) 

e  =  C^x  +  w + «  (2.40) 

z  =  C^x  +  bd^+b^^w+b^^u  (2.41) 

y  =  CyX+bdy^  +  by„w+  b^^u  (2.42) 

The  closed  loop  transfer  functions  T^a  and  for  the  mixed  problem  are  defined  as 

)*'  K  +  (2.43) 

Ted  =  Ce  )"'  B,  +  D,,  (2.44) 

The  mixed  H2/H..  controller  to  be  developed,  in  state  space  form  is  defined  as 

x^=A^x^+B,y  (2.45) 

M  =  C^x^  (2.46) 


Closing  the  loop  of  the  individual  H2  and  H„  problems  with  the  above  mixed  controller  results  in  the 
following  closed-loop  state  space  representations 
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X2  =A2Ji:2+B„w 
Z  =  C,X2+D^H' 
and 

+B^^i 
e  =  C^x^  +T)^jd 

where 


(2.47) 

(2.48) 

(2.49) 

(2.50) 

(2.51) 

(2.52) 

(2.53) 

(2.54) 

(2.55) 


r  B, 

~[BcDy,_ 

(2.56) 

C.=[C,  D^C,] 

(2.57) 

(2.58) 

Dz»-=0 

(2.59) 

^ed  -  B>ed 

(2.60) 

The  assumptions  for  the  mixed  problem  are  a  combination  of  the  assumptions  for  the  individual  H2 
and  problems.  Walker  [Wal94]  reduced  the  assumption  list  to  just  those  required  by  the  H2  problem,  so 
the  assumptions  become 
i.  D^  =  0 
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r  =  lkJL  when/^W  =  /i:2^^  .  (2.64) 

=  a  solution  to  the  H2  /  problem  for  some  7  >  7  (2.65) 

7*  =  iTed  L  -  ^mix  (2-66) 

a*s||r,J|2  when^:(^)=^:„^  (2.67) 

where  the  admissible  set  of  controllers  is  the  set  of  all  stabilizing  controllers. 


The  mixed  H2/BL  problem  can  now  be  restated  as  follows:  determine  a  K(s)  such  that 

1.  the  underlying  H2  and  IL  problems  are  stable,  i.e.,  A2  and  are  stable 


2.  y*  <y  for  some  given  y>y 

3.  IITjwll2  is  minimized. 


Walker  [Wal94]  introduced  the  following  theorem  to  further  refine  the  problem  setup. 
Theorem  1.  Let  (A^  Be  CJ  be  given  and  assume  there  exists  a  Q^=Q,l,>0  satisfying 

=0  (2.68) 

where  R  =  (7^/  )  >  0.  Then  the  following  are  equivalent: 

(i)  (A>.,Bd)  Is  stabilizable 

(ii)  A„  is  stable 
(Hi)  A2  is  stable 

Moreover,  if  the  above  hold  then  the  following  are  true: 

(iv)  lirjL<7 

(v)  the  two-norm  of  the  transfer  function  is  given  by 

lkJl2=4c,Q2C[]  =  [Q2C[cJ 

where  Q2  =  Q2  ^0  is  the  solution  to  the  Lyapunov  equation 

A2QJ  +  Q2AJ  +  B^B^  =  0 

(vi)  all  real  symmetric  solutions  (2-.  of  (2.68)  are  positive  semidefinite 

(vii)  there  exists  a  unique  minimal  solution  Q„  to  (2.68)  in  the  class  of  real  symmetric  solutions 
(viii)  is  the  minimal  solution  of  (2.68)  iff 

Re[?^  (a„  +  B^D C, )]  <  0/or  all  i 

(ix)  ||r^^||^  <(<)/ ij/9^e[/l,(A„ +B^DJj/?“‘C^ +!2„C[f?“'Cj)]<(<)  0  where  Q.„  is  the 
minimal  solution  to  (2.68). 

This  theorem  allows  the  mixed  problem  to  be  stated  as:  Determine  the  K(s)  which  minimizes  the  cost 
function 

J(A„B„C,)  =  \Tji^  =tr[Q^ClC,]  (2.69) 

where  Q2  is  the  real  symmetric,  positive  semi-definite  solution  to 
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A222  +22^2 


and  such  that 


A^Q^+Q^Ai+(Q^Cl  +B,D^,)r-'(q^CI  +B,D^,Y  +B,Bl  =0  (2.71) 

has  a  real,  symmetric,  positive  semi-definite  solution. 

The  minimization  problem  can  now  be  posed  as  a  cost  function  (2.69)  subject  to  two  equality 
constraints  (2.70)  and  (2.71).  To  solve  this  minimization  problem  a  Lagrange  multiplier  approach  is  used. 
The  Lagrangian  is 

L  =  tr[Q^C[C,]+tr^A^Q^+Q.,Al+B^Bl]x] 

(r  1  (2.72) 

+  tr{[A.2„+(2.A!;+(2„Cr  +B,BjJy| 

where  X  and  Y  are  symmetric  Lagrange  multiplier  matrices.  The  first  order  necessary  conditions  for  the 
minimum  of  this  Lagrangian  are: 

j^=X[2Qu+XiQ2+ynQat+y2Qi,=0  (2.73) 

=  XJ^Q.Cl  +  X^Qj^Cl  +  €]_  +  Y^Q^Cl 

+  +y2BcV>,  +(n22.  +r2GjJcJ'M+(T,^e„,  +Y^Q,)clDiM  =0  (2.74) 

■^‘*2i2  +  2^1222  +  ^I22|2  +2?2CcG2  +  B^_Y^2Qb 

+  R^b  Qa  yi  Qab  +  RL  Qa  XuQb R^b Qab ^^12 Qab  +  RL  Qab ^2  Qb 

+  Rb^cQaby\Qab  Rb^cOby^Qab  Rb^cQabyuQb  Rb^cQbyiQb 

+Px{Y,Q,b+yuQb)+P2iyuQab+y2Qb)=^  (2.75) 


■^  =  A2G2+22A2  +BwBw  =0 


—  =  A[X+XA2+C[C,=0 


^  =  A„G„+G„Ai+(G^C^+B^Drd)^“‘(G„C[+B^D[j^+B^Bj  =0 
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^=(a„  +b,dL/?-'c,  +q„cIr-^cJ  y 
+  y(a^  +BXd^“'c,  +2^cf^-'cJ=o 


where 


m  =  r-'d,,dI, 


P,=DlR-^D,,Bl 


P2=DlMBl 


Q^  = 


Q\  Qn 
Qii  Qi 


x= 


A,  X 

■T 
12 


X,^,  X 


= 


Qa  Qab 

Qab  Qb 


Y  = 


n  Yn 
Yn  J^2, 


b„b;,= 


^c^yw 


[bI  dI,bI 


n  YuBl 
BcVil  b.v^bJ 


B,(DLi?-'D,,  +/)bJ  =1  +/)[5j 


BcV,B^ 


clc,  = 


.^c  ^zuj 


[Cz 


^12  Q 


c^r!^  cJr^c, 


(2.79) 

(2.80) 
(2.81) 
(2.82) 

(2.83) 

(2.84) 

(2.85) 

(2.86) 

(2.87) 


(2.88) 


(2.89) 
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(2.90) 


_c  ^eu  _ 

_  ^ab^c 

~  C^Rl  CjR,C, 

Equations  (2.76)  and  (2.78)  are  the  original  constraint  functions,  equations  (2.70)  and  (2.71).  There 
are  two  possible  solutions  to  (2.79).  One  solution  to  (2.79)  is  that  Y  =0;  the  other  possibility  is  that 

(a„  +2ooCjf?~’Cj)  is  neutrally  stable.  The  first  solution,  Y  =0,  corresponds  to  the  an 

inactive  He  constraint.  The  second  solution  implies  that  solution  lies  on  the  boundary  of  the  He  constraint 
and  Qe  is  the  neutrally  stabilizing  solution  for  the  BL  Riccati  equation  (2.71).  When  Y  =0,  the  Lagrangian 
reduces  to 

i.  =  fr[(22C[cJ+rr{[A2e2+!22A[+B„B;:]x}  (2.91) 

which  is  the  same  as  the  H2  Lagrangian.  Therefore,  for  Y  =0,  Kn,ix  =  K2opf  Walker  [Wal94]  developed  the 
following  theorem: 

Theorem  2  Assume  tic  ^  where  tic  is  the  order  of  the  controller  and  nz  is  the  order  of  the  plant  for  the 
H2  problem.  Then  the  following  hold: 

1.  Ify  <%no  solution  to  the  mixed  H2/H„  problem  exists 

2.  If  y  <j  <J ,  is  such  that  7^=7 

3.  Ify>y,K2^  is  the  solution  to  the  mixed  H2/H„  problem. 

When  the  order  of  the  controller  is  chosen  greater  than  or  equal  to  that  of  the  H2  problem,  the 
solution  to  the  mixed  H2/H00  problem  will  lie  on  the  EL  constraint  boundary.  Theorem  2  implies  that  each 
point  on  the  mixed  H2/EL  curve  occurs  on  the  H„  constraint  boundary',  y*=T>  for  y  <  y  <  y .  In  this  region. 


a*  is  a  monotonically  decreasing  function  of  y.  A  typical  a  versus  y  curve  is  shown  in  Figure  2.7. 


Figure  2.7.  Typical  Mixed  H2/H„  a  Versus  y  Curve 

2.2.  Handling  Qualities  Analysis 

There  are  several  handling  qualities  prediction  criteria  commonly  used  to  evaluate  flight  control 
systems.  These  criteria  are  based  on  historical  data  of  pilot  ratings,  and  are  outlined  in  the  Department  of 
Defense  MIL-STD-1797A,  Flying  Qualities  of  Piloted  Aircraft  [Mil90].  Two  criteria  were  used  to  evaluate 
the  flight  control  systems  designed  by  the  various  methods  in  this  thesis.  The  primary  criteria  used  was 
Hoh’s  bandwidth  criteria.  The  TPS  project  this  thesis  was  based  on,  uses  this  criteria.  The  Ralph  Smith 
(RSmith)  criteria  was  also  considered,  and  is  the  criteria  TPS  is  currently  using  on  an  updated  flight  control 
system  project  [TPS95].  The  next  two  subsections  briefly  describe  these  handling  qualities  criteria.  For  a 
full  description  see  the  [MIL90],  [Hoh96]  and  [SG78]. 

2.2.1.  Bandwidth  Criteria 

The  bandwidth  requirement  is  based  on  the  hypothesis  that  if  the  aircraft  has  good  response 
characteristics  over  a  sufficiently  wide  range  of  pilot  control  input  frequencies,  then  the  aircraft  will  have 
favorable  handling  qualities  [HohSl].  This  evaluation  method  examines  two  parameters  -  an  equivalent 
time  delay,  to  account  for  the  higher  order  dynamics  of  the  aircraft,  and  the  bandwidth  of  the  pitch  attitude 
to  stick  force  (or  displacement)  transfer  function. 

The  bandwidth  is  defined  as  the  lesser  of  two  frequencies,  and  .  The  bandwidth 

due  to  phase  margin,  co  aWp^ase  ’  defined  as  the  frequency  at  a  phase  margin  of  45  degrees.  The  frequency 
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where  the  phase  angle  is  -180°,  labeled  ooiso,  and  the  gain  which  corresponds  to  that  frequency  are 
determined.  The  bandwidth  due  to  gain  margin,  ^  frequency  at  the  gain  which  is  6 

dB  above  the  gain  amplitude  corresponding  to  coigo-  The  bandwidth  determination  is  graphically  depicted  in 
Figure  2.8. 


Figure  2.8.  Definition  of  Bandwidth  Frequency  (Obw  from  Open  Loop  Frequency  Response 


The  second  parameter  needed  to  determine  predicted  handling  qualities  level  is  the  equivalent  time 
delay,  Xp.  This  parameter  is  computed  using  the  formula 


(*l’2a)i8o + 

2(57.3)(co,go) 


(2.92) 


where  COigo  is  the  frequency  corresponding  to  -180°  phase  and  (t>2Q)]go  phase  angle  at  a  frequency 

twice  as  great  as  the  frequency  where  -180  degrees  of  phase  angle  occurs. 

These  two  parameters  are  plotted  on  the  CObw  versus  Xp  Bandwidth  requirements  diagram  to 
determine  the  predicted  handling  qualities  level.  There  are  separate  diagrams  for  different  flight  phases. 
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(Obw  (rad/sec) 

Figure  2.9.  Bandwidth  Requirements,  Category  C  Flight  Phase 

Figure  2.9  depicts  the  level  boundaries  for  the  bandwidth  criteria,  for  Category  C  Flight  Phase.  This  thesis 
is  only  concerned  with  approach  and  landing,  Category  C  Flight  Phase,  so  only  this  diagram  is  reproduced 
here. 

2.2.2.  Ralph  Smith  Criteria 

A  pitch  and  roll  axis  handling  qualities  criteria  was  proposed  by  Ralph  Smith  and  Norman  Geddes 
[SG78].  The  RSmith  criteria  predicts  a  pilot-vehicle  cross-over  frequency,  a  Cooper-Harper  rating,  and  a 
pilot-induced  oscillation  (PIO)  rating  in  both  the  pitch  and  roll  axes.  Only  the  pitch  axis  portion  of  this 
criteria  will  be  presented  here.  The  proposed  pitch  axis  requirements  are  restricted  to  longitudinal  mode, 
short  period  dynamics  of  Class  FV  aircraft  (high-maneuverability  -  fighter  type),  Flight  Phase  Category  A 
(non-terminal  flight  phases  that  require  rapid  maneuvering,  precision  tracking,  precise  flight-path  control). 
The  requirements  proposed  in  [SG78]  were  derived  from  flight  test  data  for  a  variety  of  flight  control  tasks; 
the  power  approach  task,  however,  was  not  included  in  this  data  base.  The  Smith-Geddes  paper  does  argue 
that  this  criteria  is  “sufficient  to  ensure  acceptability  of  front-side  power  approach  handling  qualities . . .  [but] 
may  be  too  stringent  for  direct  use  as  design  specifications  for  the  power  approach  condition.”  [SG78:142] 
The  test  bed  used  for  this  research  is  a  Class  II  aircraft,  but  it  is  simulating  Class  IV  aircraft  dynamics. 

Also,  since  TPS  is  currently  using  the  RSmith  criteria  to  evaluate  the  approach  and  landing  task  in  the 
current  TPS  FCS  project,  the  RSmith  criteria  was  considered  during  design  as  a  secondary  prediction  of 
handling  qualities. 
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The  RSmith  criteria  for  the  pitch  axis  can  be  summarized  in  a  six  step  process  [TPS95].  It  requires 
measuring  a  slope,  determining  a  criterion  frequency,  and  obtaining  two  phase  angles.  Each  of  these 
parameters  has  an  associated  handling  qualities  requirement.  The  six  steps  and  the  associated  requirements 
are  as  follows: 

Step  1:  From  a  pitch  rate  time  history,  determine  the  elapsed  time  to  the  first  peak,  tg,  resulting  from  a  step 
input.  The  time  to  first  peak  must  be: 

Level  1:  0.2<  r,  <  0.9 
Level  2:  none  proposed 


Step  2: 


Determine  the  slope  parameter,  S,  the  slope  of  the  Bode  magnitude  curve  of 


dB/octave,  between  1  and  6  rad/sec.  In  equation  form,  S  is  defined  as 


S  =  - 


dca 


-U<^c) 


dB  /  octave 


(2.93) 


The  handling  qualities  requirement  is 

Level  I :  S  <  -2  dB/octave 
Level  2:  none  proposed 


Step  3;  Determine  the  pilot- vehicle  cross-over  frequency,  also  called  the  criterion  frequency,  COb,  using  the 
formula: 


(Ob  =  6.0  +  0.24  5  (radians/second) 


(2.94) 


where  S  is  the  slope  of  the  Bode  magnitude  curve  of 


,  determined  in  step  2. 


Step  4:  Determine  the  phase  angle  of  - (Mc)  use  Figure  2. 1 0  to  determine  the  predicted 


Cooper-Harper  rating. 
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Step  5:  If  Z-— (jco^)  >  -180°,  susceptibility  to  pitch  attitude  PIO  is  not  predicted.  However,  if 


O  /  V 

Z - (y®c  j  -  -180°,  a  PIO  may  result  from  pitch  attitude  tracking. 


Figure  2.10.  Criterion  Phase  Angle  Versus.  Cooper-Harper  Rating 


Step  6:  Determine  the  normal  acceleration  phase  angle  parameter  ‘I>(j(0b)  from  the  following  formula: 

®(ycoJ  =  Z-^{y®J-14.3co,  (2.95) 

*  es 

where  is  normal  acceleration  at  the  pilot’s  location  in  the  airplane.  The  -M.SCOb  term  is  the 

phase  lag  introduced  by  typical  pilot  dynamics  and  is  modeled  by  a  transport  delay  of  0.25 
"z 

seconds.  The  angle  of -^(yco^)  is  measured  from  a  Bode  phase  angle  curve.  For  conventional 

^es 

flight  control  systems,  the  PIO  rating  can  be  predicted  using  <I>(j(iOb)  in  Figure  2.1 1. 
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6 


Level  3 
or 

W  orse 

Level  2 

1  1 

T - , - , - j - 1  I - 

150  -170  -190  -210  -230  -250  -270 

<P  {j(i)c),  degrees 

Figure  2.1 1.  Predicted  Pitch  Axis  PIO  Rating  for  Conventional  Flight  Control  Systems 

These  requirements  should  be  considered  a  set  of  necessary  conditions.  Each  must  be  satisfied  at  the 
Level  1  or  2  values  for  the  overall  handling  qualities  rating  to  be  Level  1  or  2,  respectively.  If  any  one  is 
violated,  the  overall  rating  is  equal  to  the  level  for  which  the  violation  exists. 

2.3.  TPS  Project 

The  US  Air  Force  Test  Pilot  School  (TPS),  as  part  of  its  curriculum,  requires  the  students  to 
complete  a  Flight  Control  System  (FCS)  design  project.  This  project  includes  a  scenario  involving  analysis 
of  a  proposed  contractor  designed  flight  control  system,  then  a  redesign  of  the  FCS  based  on  two  specified 
design  structure  options.  The  student  test  teams  choose  one  of  two  design  options.  The  structure  of  these 
design  options  (i.e.  the  designated  feedback  loops)  is  fixed,  and  the  test  team  is  limited  to  selecting  FCS 
gains  and  stick  characteristics.  The  FCS  is  for  a  Class  IV  type  aircraft  performing  the  approach  and  landing 
task.  The  project  specifies  evaluating  the  handling  qualities  of  the  design  by  applying  Hoh’s  Bandwidth 
Criteria.  The  FCS  design  is  then  implemented  in  the  Calspan  Variable  Stability  Simulator  (VSS)  Learjet 
and  the  handling  qualities  are  evaluated,  using  the  Cooper-Harper  rating  scale,  by  the  student  test  pilots 
performing  a  lateral  offset  approach  to  a  spot  landing.  This  project  was  used  as  a  foundation  for  this  thesis 
research.  Only  the  portion  of  the  TPS  FCS  project  used  for  this  research  will  be  presented  here.  Design 
option  one  was  chosen  as  the  basic  structure  for  the  flight  control  systems  designed  in  this  thesis.  The 
system  was  restricted  further  by  fixing  the  stick  dynamics  used  to  the  baseline  stick  dynamics  and  feel 
system  presented  in  the  TPS  project. 
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The  TPS  PCS  project  strictly  uses  classical  design  methods  (i.e.  Root  Locus,  Bode  frequency,  and 
time  response  analysis)  and  transfer  function  models  of  the  aircraft  dynamics.  The  pitch  angle  to  elevator 
deflection  and  the  angle  of  attack  to  elevator  deflection  transfer  functions  specified  in  the  TPS  project  are: 


0(5)  2.85(j+0.06)(5  +  0.87) 

5(^  “  (5  + 155)(^ - 0.38)(5 + 0.06 ±  0.23 ;) 

a(j)  0.04l(j+0.016±0.21  j){s+10.2) 
+ If  5)(5 - 0.38)(5 + 0.06 ±  0.23 j) 


(2.96) 


(2.97) 


The  structure  of  design  Option  1,  shown  in  Figure  2.12,  is  based  on  an  inner  angle  of  attack  feedback  loop 


(with  feedback  gain  of  KA,  and  a  low  pass  noise  filter),  an  outer  pitch  rate  feedback  loop  (with  feedback 


gain  KQ),  and  a  forward  command  gain  (KFL).  A  nominal  value  and  a  range  for  these  gains  was  also 


specified,  and  are  shown  in  Table  2.1. 


DISPLACEMENT/FORCE 


FEEL  SYSTEM 


BREAKOUT 

GRADIENT 

F  — » 

D/^_l  A 

DFG=0.17 

es  ^ 

dUsI.U 

r* 

(IN) 

(lbs) 

16^ 


i^+2*a7*16s+I6^ 


COMMAND 

GAIN 


FORWARD  LOOP 
COMPENSATOION 


ACTUATOR 


(IN) 


KCal.O 


(DEG)' 


0<KC^O 


KFULO 


(DEG) 


0<KFL^ 


20 

^  +  20 


5. 

(DEG) 


q-GAIN 


Q 

(DEG/SEC) 


a 

(DEG) 


J  KQ=:L0 


0<;kqs2.5 
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FILTER 
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j+10 

AOA- 
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KA=2^ 

(DEG) 
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Figure  2.12.  TPS  FCS  Design  Option  1  (Nominal) 


GAIN 

NOMINAL  VALUE 

RANGE 

KFL 

1.0 

0<KFL<5 

KA 

2.5 

0<KA<5 

KQ 

1.0 

0  <  KQ  <  2.5 

Table  2.1.  TPS  FCS  Design  Option  1,  Nominal  Gains  and  Ranges 
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The  general  structure  of  this  TPS  PCS  project  was  the  foundation  for  the  flight  control  systems 
designed  in  this  thesis.  The  classically  designed  system  followed  this  setup  exactly.  The  modern  designs 
(H2,  KL,  and  mixed  H2/H„  designs)  only  fed  back  pitch  rate  and  angle  of  attack,  like  the  TPS  project,  but 
used  state  space  methods  and  did  not  limit  the  controller  to  simple  gains  and  a  low  pass  filter.  The 
controller  was  also  not  limited  to  negative  feedback,  nor  to  the  feedback  path. 


2.4.  CALSPAN  Variable  Stability  Simulator  II  Learjet 

This  thesis  research  implemented  the  PCS  designs  on  the  CALSPAN  Variable  Stability  Simulator  II 
(VSS  n)  Learjet,  tail  number  N102VS.  The  VSS  II  is  a  modified  Learjet  Model  25  that  functions  as  a  three 
axis  in-flight  simulator  [Bal93].  It  is  operated  by  CALSPAN  under  an  experimental  license  from  the  PAA. 
The  cockpit  has  a  set  of  side-by-side  controls.  The  control  yoke  at  the  left  seat,  for  the  safety  pilot, 
maintains  the  Learjet’ s  conventional  flight  control  system.  The  evaluation  pilot’s  center  control  stick  in  the 
right  seat  uses  a  fly-by-wire  response  feedback  system.  The  variable  stability  and  variable  control  system 
consist  of;  variable  feel  system,  aircraft  motion  sensors  and  associated  signal  conditioning,  control  system 
simulation  computer,  control  surface  servos,  digital  configuration  control  system,  engage/disengage  and 
safety  monitor  logic,  and  recording/playback  capability.  The  system  can  be  disengaged  by  any  of  four 
manual  disconnect  switches  which  are  installed  on  each  of  the  three  control  sticks/yokes  and  on  the  glare 
shield.  In  addition,  automatic  safety  trips  are  provided.  In  the  event  of  incapacitation  of  the  safety  pilot  or 
certain  control  cable  failures,  the  aircraft  can  be  flown  by  the  evaluation  pilot  as  a  nearly  normal  Learjet 
using  the  VSS  in  the  “Ply-by-Wire”  (PBW)  mode.  All  basic  Learjet  systems  are  available  in  the  PBW 
mode  except  for  nose  gear  steering.  The  handling  qualities  are  those  of  the  basic  aircraft  with  the  yaw 
damper  on.  There  are  no  safety  trips  in  the  emergency  PBW  mode. 

Hydraulic  power  for  the  variable  stability  actuators  is  obtained  from  the  existing  Learjet  hydraulic 
system  which  provides  four  gallons  per  minute  (gpm)  per  engine.  Estimated  maximum  flow  demand  to 
operate  all  servo  actuators  is  3.35  gpm.  Maximum  demand  for  normal  Learjet  flaps,  spoilers,  gear  and 
brakes  is  under  four  gpm.  Solenoid  operated  valves  to  the  variable  stability  actuators  are  designed  fail  safe 
to  prevent  hydraulic  locks  on  the  actuators. 
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Some  operating  limitations  that  apply  to  the  Leaqet  are  listed  in  Table  2.2. 


VSS  OFF 

VSS  ON 

SPEED 

LIMITATIONS 

356  KIAS 

0.82  Mach 

325  KIAS 

G  LIMITATIONS 

+4.4  to  -1.0 

+2.8  to  +0.15 

Table  2.2.  Learjet  Limitations 


The  control  system  simulation  computer  uses  MATLAB®  SIMULINK  to  modify  the  perceived 
aircraft  dynamics.  For  the  longitudinal  axis,  only  short-period  dynamics  are  simulated.  While  the  system 
programming  and  inter-workings  are  proprietary  to  Calspan,  the  following  guidance  for  implementing  FCSs 
was  provided: 

1.  Suggest  use  of  MATLAB®  version  4.2C.1  and  SIMULINK  1.3  or  later. 

2.  Recommend  controllers  be  cast  in  pole/zero/gain  form  using  scalar  parameters  and 
integrators.  Do  not  use  built  in  transfer  function  blocks  unless  the  parameters  will  never 
change.  An  example  block  diagram  is  shown  in  Figure  2.13. 

3.  Do  not  use  MATLAB®  function  blocks. 

4.  Setup  file  (MATLAB®  script)  must  contain  only  constants,  no  expressions. 

5.  Do  not  use  masked  blocks  with  dummy  parameters. 

6.  Do  not  use  a'^b,  use  pow(a,b) 

7.  Do  not  use  matrix  SIMULINK  parameters  unless  they  won't  ever  change 

8.  Use  named  scalar  parameters  for  anything  you  might  change  and  initialize  it  with  your 
setup  file. 

9.  Don't  forget  initialization  of  integrators.  A  "system  engage"  discrete  is  available  for 
controlling  reset  integrators. 

10.  Use  Euler  integration,  with  a  0.01  step  size. 

11.  Please  make  sure  your  top  level  interface  is  well  defined  including  engineering  units, 
etc.  Inputs  are  generally  pilot  forces  or  stick  displacements  and  engage  status  logicals 
and  sensor  signals.  Outputs  are  generally  just  surface  actuator  commands. 


The  transfer  function. 


^  ec  ^  loop 


^  (t„s+z) 

K.5„  — Ua 


(XjS+p) 
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3)  t-K±l  - 
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Hii— I 


Figure  2.13.  Example  SIMULINK  Block  Diagram  for  a  VSS  n  Acceptable  System 


From  the  list  of  guidelines  provided  by  Calspan,  the  critical  factors  were  that  the  system  must  be 
stable  using  Euler  integration  with  a  0.01  seconds  step  size,  and  the  integrators  must  by  able  to  be  initialized 
and  zero-ized  (reset)  in-flight  with  the  system  running  -  therefore  the  SIMULINK  state-space  blocks  were 
not  directly  implementable  in  the  VSS  II  system.  How  this  was  accomplished  will  be  covered  in  Chapter  4, 
Implementation. 

2.5.  Summary 

This  chapter  has  covered  the  background  theory  used  in  this  thesis  and  background  information 
concerning  the  resources  and  the  foundation  design  set-up  on  which  this  thesis  is  based.  This  chapter 
started  with  classical  design  theory  fundamentals,  and  then  described  the  basics  of  Hg  and  IL  optimization. 
Next,  mixed  Hg/HL.  optimization  theory  was  briefly  discussed.  Two  handling  qualities  design  criteria  were 
explained;  specifically  addressed  were  Hoh’s  Bandwidth  criteria  and  RSmith  criteria.  The  TPS  project, 
which  was  the  foundation  for  this  research,  was  introduced.  Finally  the  Calspan  VSS  II  Learjet  was 
described,  including  the  design  implementation  guidelines  specified  by  Calspan.  The  next  chapter  will 
address  the  set-up  used  to  formulate  the  FCS  designs. 
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III.  Model  Set-up  and  Design 


3. 1.  Aircraft  Plant  Model 

The  TPS  PCS  project  provided  an  aircraft  plant  model  in  transfer  function  form.  This  model  only 

(  0(^)  "1  (  «(^)  1 

provided  the  pitch  angle  — ^  and  angle  of  attack  — — -  dynamics  (see  equations  2.96  and  2.97).  The 

modem  control  theory  designs  investigated  in  this  thesis  required  a  state-space  description  of  the  aircraft 
dynamics  [RB86].  Calspan  provided  a  set  of  stability  derivatives  and  flight  conditions  for  the  aircraft 
model  used  in  the  TPS  PCS  project  during  the  past  few  years.  This  list  of  stability  derivatives  and  flight 
conditions  and  the  MATLAB™  script  file  used  to  calculate  the  state-space  model  are  shown  in  Appendix  A. 
At  the  time  the  state-space  model  was  created,  the  VSS  11  limitation  of  only  simulating  short  period 
dynamics  was  not  known.  Therefore,  since  it  is  more  complete,  a  full  state-space  model  was  determined 

based  on  the  data  provided  by  Calspan.  The  states  used  were:  q  -  pitch  rate  (degrees/second),  CL  -  angle  of 

attack  (degrees),  v  -  forward  true  velocity  (feet/second),  y-  glide-path  angle  (degrees),  and  once  the 

actuator  was  included  in  the  model,  x^ct  -  the  actuator  state.  The  sign  convention  for  elevator  deflection 
used  by  TPS  is  the  NASA  sign  convention,  a  positive  elevator  deflection  implies  the  trailing  edge  moves 
down  and  will  produce  a  negative  pitching  motion  [TPS95].  The  actuator  model  used  in  the  TPS  project, 

20 

- ,  was  added  to  the  aircraft  plant  model  for  the  design  and  evaluation  process.  The  variables 

S-H20 

considered  for  analysis,  the  output  variables,  were  the  states,  CC,  q,  J,  v,  as  well  as  6  -  pitch  angle  (degrees), 
and  nj  -  normal  acceleration  (g’s).  The  resulting  state-space  model  was 


i:  =  Ax-fBM 

(3.1) 

y  —  Cx  +  Du 

(3.2) 

where 

'q 

a 

x=  V 
Y 

f'act 
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«  =  [S.c] 


-0.38 

0J8 

0.00 

0.00 

2.86 

1.00 
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0.00 

0.04 

0.00 
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-056 

0.00 

0.00 

0.85 
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0.00 
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0.00 

0.00 

0.00 

0.00 

-20.00 
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0.00 

0.00 
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0.00 
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1.000 

0.000 

0.000 

0.000' 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1.000 

0.000 
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0.000 
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0.000 
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0.000 

0.000 

0.000 

1.000 

0.000 

0.000 

0.000 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1.000 

D  = 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


This  state-space  model  can  be  transformed  back  to  transfer  function  form  to  verify  accuracy  of  the 
model  (compared  to  the  transfer  functions  given  in  the  TPS  project)  and  for  use  in  the  classical  design  as 
well  as  evaluation  of  the  designs.  The  transfer  function  representation  of  this  state-space  model  is 
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0(^)  2.86(20){5+0.06)(j+0.84) 

^  ( j + 20)(5 + U3)(5  -  0.39)( J + 0.07  ±  0.23 ;) 

a(s)  0.04  l(20)(j+ 0.024  ±0.21 7)(^+70.2) 

^  ~  {s+ 20)(j  +  1.53)(5 - 0.39)( j + 0.07  ±  0.23 ;) 

When  compared  to  (2.96)  and  (2.97),  the  differences  are  small. 


3.2.  Wind  and  Sensor  Noise  Model  /  Intensity  Determination 

One  of  the  objectives  of  the  flight  control  designs  is  to  provide  stability  and  performance  in  the  face 
of  noises  acting  on  the  system.  Both  for  design  and  analysis,  the  noises  needed  to  be  categorized  and 
modeled.  The  two  types  of  noise  considered  were  wind  gust  noise  and  sensor  noise.  The  limited  nature  of 
the  flight  test  prevented  the  noise  rejection  characteristics  of  the  designs  to  be  quantified  during  the  flight 
test.  Prior  to  flight  test,  computer  analysis  was  performed  to  evaluate  each  design’s  noise  rejection 
characteristics. 

The  wind  gust  noise  can  be  characterized  as  a  random  low  frequency  disturbance  in  angle  of  attack, 
a.  It  is  modeled  as  an  angle  of  attack  perturbation  by  inserting  the  deviation  directly  into  the  core  plant 
(bypassing  the  actuator),  see  Figure  3.1.  This  is  done  by  multiplying  the  filtered  wind  times  the  second 
column  of  the  A  matrix  of  the  aircraft  plant  and  adding  this  to  the  states.  The  second  column  of  A  will  be 
denoted  by  T.  For  the  purposes  of  design,  the  wind  noise  is  modeled  as  a  zero-mean.  Unit  intensity  White 
Gaussian  Noise  (UWGN),  with  no  filtering  or  ‘coloring’  added  since  the  filtering  proves  only  to  increase 
the  order  of  the  controller  and  not  improve  the  controller’s  performance.  The  design  process  included 
weighting  the  noise  input  to  obtain  the  desired  noise  rejection  results.  For  evaluation  of  the  designs,  the 
UWGN  with  a  sample  time  of  0.01  seconds,  was  passed  through  a  low'  pass  filter,  represented  by 


W„  =■ 


0.1 


(3.5) 


*  ^-(-0.5 

This  filter  removes  the  high  frequency  component  of  the  UWGN  to  result  in  a  ‘wind  like’  disturbance  of 
reasonable  intensity  for  approach  and  landing.  However,  when  comparing  designs,  to  accentuate  the 
differences  in  designs  caused  by  the  ‘wind  noise’,  the  intensity  was  increases  by  a  multiple  of  five.  This 
proved  beneficial  for  analysis. 
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Sensor  noise  can  be  characterized  as  a  high  frequency  disturbance  entering  the  system  at  the  output 
of  the  plant,  as  a  result  of  measuring  the  signals.  The  measurements  are  corrupted  by  UWGN  with  a 
strength  of  0.1  deg^-sec  and  a  sample  time  of  0.01  seconds.  For  design,  the  strength  of  these  noises  were 
nominally  set  at:  =  0.01, for  the  noise  intensity  of  the  alpha  measurement,  and  Nq=0.001,  for  the  pitch 
rate  measurement.  The  rate  measurement  was  expected  to  be  better  then  the  angle  measurement,  but  no 
actual  truth  data  on  the  system  to  be  tested  was  available.  These  values  were  then  tuned  with  the  output 
weightings  to  obtain  the  desired  design.  The  evaluation  model  used  a  high  pass  filter  to  shape  the  UWGN 
input.  The  filter  used  was 


W.  = 


0.1(5  +  10) 


5  +  100 

with  a  sample  time  of  0.01  seconds.  Figure  3.1  depicts  the  evaluation  full  model,  including  the  noise 
models  used  to  evaluate  the  designs  for  noise  rejection. 
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Figure  3.1.  General  Evaluation  System 


3.3.  Evaluation  Models 

Each  of  the  controllers  designed  were  evaluated  by  considering  the  time  response  generated  by 
SMULINK  simulations,  handling  qualities  predictions  based  on  the  Bandwidth  criteria  (RSmith  criteria 
was  also  considered),  singular  value  frequency  plots,  and  vector  stability  margins.  The  results  from  these 
evaluations  are  discussed  in  the  individual  design  models  and  resulting  design  sections  that  follow,  with  the 
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time  history  plots  presented  in  Appendix  C.  The  latter  two  considerations  were  used  to  expedite  the  design 
process,  and  were  accomplished  by  a  simple  series  of  MATLAB™  commands.  The  MATLAB™  script  fdes 
used  are  presented  in  Appendix  B.  The  next  two  sub-sections  present  the  set-up  for  the  time  history  and  the 
handling  qualities  analysis. 

3.3.1.  Time  Domain  Evaluation 

The  MATLAB™  SIMULINK  simulation  program  was  used  to  obtain  time  histories  of  the  feedback 
control  systems  designed.  The  same  SIMULINK  diagram,  shown  in  Figure  3.1,  was  used  for  all  the  modern 
design  methods  to  obtain  the  simulated  time  histories.  The  classical  design  used  the  SIMULINK  diagram 
shown  in  Figure  3.2,  to  obtain  the  simulated  time  histories.  Different  diagrams  were  used  to  accommodate 
the  location  of  the  controller  in  the  feedback  path. 


Figure  3.2.  Classical  Evaluation  Model 


All  the  SIMULINK  simulations  used  the  LINSIM  numerical  integration  routine,  a  time  step  of  0.01 
seconds,  and  for  a  duration  of  10  seconds.  Both  a  step  input  and  an  initial  condition  of  5°  angle  of  attack 
were  considered,  with  no-noise  introduced  and  with  five  times  the  nominal  noise  intensity.  The  step  input 
was  used  to  evaluate  the  tracking  response  of  the  system  and  the  initial  condition  was  used  to  evaluate  the 
perturbation  rejection  characteristics  of  the  system.  Both  input  types  were  used  to  consider  the  noise 
rejection  properties  of  the  designs. 
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3.3.2.  Handling  Qualities  Evaluation 

The  handling  qualities  rating  of  the  flight  control  systems  designed  were  predicted  using  Hoh’s 
Bandwidth  criteria  as  primary  evaluation  method  and  RSmith’s  criteria  as  a  secondary  evaluation  method. 
The  Interactive  Flying  Qualities  Toolbox  for  Matlab,  version  0.01  Beta,  February  1995  was  used  to  predict 
the  handling  qualities  of  the  systems.  This  toolbox  takes  as  its  input  the  pitch  angle  to  stick  force  transfer 
function.  The  nominal  stick  dynamics  defined  in  the  TPS  FCS  project  were  used  to  define  this  closed  loop 
transfer  function. 

3.4.  Design  Models  and  Resulting  Designs 
3.4.1.  Classical  Design  Set-up  and  Resulting  Design 

The  classical  design  was  included  in  this  thesis  as  a  control  case  for  the  flight  test  and  as  a  means  to 
validate  the  bare  airframe  aerodynamics  simulated  by  the  Calspan  VSS  n  Learjet.  The  set-up  and  structure 
was  based  on  Design  Option  One  of  the  TPS  FCS  project  described  in  Chapter  II.  The  system  can  be 
represented  by  the  block  diagram  shown  in  Figure  3.3.  The  state-space  aircraft  dynamics  model  shown  in 
Figure  3.3  incorporates  the  actuator  specified  by  the  TPS  project.  It  used  negative  feedback  of  angle  of 
attack  and  pitch  rate,  with  the  controller  gains  in  the  feedback  path  and  the  alpha  channel  incorporating  a 
low  pass  filter.  The  ranges  on  the  three  variable  gains  were  set  by  the  TPS  project. 


This  made  a  vary  iterative  process  manageable  to  an  inexperienced  designer.  The  nice  aspect  of  the 
classical  design  method  with  this  simple  structure  is  the  physical  interpretation  of  the  variables  and 
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feedback  loop  gains.  Increasing  will  generally  increase  the  natural  frequency  and  decrease  the  damping 
ratio  of  the  second-order  aircraft  response  modes.  Increasing  Kq  generally  will  increase  the  damping  ratio  of 
the  system  [Bal93].  Determining  the  balance  of  the  two  which  meets  the  desired  handling  qualities  and 
performance  constraints  is  the  challenge  to  the  designer.  There  is  no  way  to  determine  if  the  gains  chosen 
provide  either  the  highest  stability  margins  possible,  or  the  best  noise  rejection  characteristics  obtainable, 
for  a  given  performance  objective  and  handling  qualities  rating.  The  controller  gains  or  states  do  maintain 
their  physical  interpretation  throughout  the  process.  The  gains  chosen  were 

Ka  =  3.4 

Kq  =  2.3 

Kfl  =  1.0 

or  written  in  state-space  form,  the  controller  can  be  represented  by 

Ak  =  [-10.0]  Bk  =  [1.0  0.0] 

Ck  =  [34.0]  Dk  =  [0.0  2.3] 

The  time  histories  are  shown  in  Appendix  C.  The  handling  qualities  analysis  results  for  the  Bandwidth 
criteria  are  shown  in  Figure  3.4.  and  the  results  for  the  RSmith  criteria  are  shown  in  Table  3.1. 


Bandwidth  (rad/sec.) 

Figure  3.4.  Category  C  Bandwidth  Requirements:  Classical  Design 
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Parameter 

Value 

Level 

S  (dB/octave) 

-3.97 

1 

5.05 

-150.84 

2 

0.48 

1 

Table  3.1.  RSmith  Handling  Qualities  Predicted  Level  for  the  Classical  Design 


The  structure  of  this  system  contains  two  feedback  loops,  Ot  and  q,  but  only  one  input,  5^.  To  fully 

describe  the  stability  margins,  four  types  of  vector  stability  margins  were  calculated.  The  first  three  are 
single-input-single-output  (SISO)  methods,  considering  perturbations  in  one-loop-at-a-time,  closing  the 
remaining  feedback  loop.  Input  vector  margins  and  output  vector  margins  were  determined  in  this  way.  It 
should  be  emphasized  that  vector  margins  are  conservative  bounds  on  the  allowable  system  perturbations, 
and  that  conventional  stability  margins  are  guaranteed  to  be  as  good  or  better.  The  fourth  type  of  vector 
stability  margin  determined  was  a  MIMO  vector  margin.  It  considers  breaking  both  loops  simultaneously 
and  is  the  most  conservative  bound  on  the  vector  gain  and  phase  margins.  For  a  more  complete  discussion 
of  vector  stability  margins  in  general,  see  [FPE94:422-425]  and  for  a  derivation  of  MIMO  vector  margin 

calculations  see  [RB86].  Table  3.2  presents  the  vector  stability  margins  for  the  classical  design.  It  should 

\ 

be  noted  that  the  TPS  FCS  project  did  not  address  or  specify  stability 


Margin  Type 

Vector  Gain  Margin 
(dB) 

Vector  Phase  Margin 
(degrees) 

Input  Margin 

[5.86  -oo] 

±  60.84 

Output  Margin 
(Ot  loop  ‘broken’) 

[  11.70  -8.03] 

±43.43 

Output  Margin 
(q  loop  ‘broken’) 

[  6.87  -8.66  ] 

±  36.78 

MIMO  Margins 
(sensitivity) 

[  3.03  -2.24  ] 

±  17.93 

MIMO  Margins 
(complementary  sensitivity) 

[  2.36  -3.24  ] 

Table  3.2.  Classical  Design  Vector  Stability  Margins 


margins  of  any  kind.  Therefore  there  were  no  minimum  acceptable  margin  levels  targeted.  The  MIL-STD 
[Mil90]  refers  to  conventional  stability  margins  of  ±6  dB  gain  margin  and  ±45°  phase  margin.  Since  vector 
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margins  are  conservative,  values  approaching  those  in  the  MIL-STD  were  considered  safe  to  flight  test. 

The  MEMO  vector  margins  were  considered  only  for  comparison  between  designs. 

3.4.2.  H2  Design  Set-up  and  Resulting  Design 

The  H2  design  problem  was  set-up  as  a  one  degree  of  freedom,  angle  of  attack  command  controller  in 
the  feed-forward  loop  with  output  weightings  to  influence  the  resulting  controller  and  a  weight  on  the  input 
to  penalize  the  control  power  of  the  system.  Figure  3.5  depicts  the  block  diagram  of  the  H2  problem 


Figure  3.5.  H2  Design  Model  Block  Diagram 


set-up.  Here  the  actuator  is  separate  from  the  aircraft  plant,  and  the  matrix  H  represents  the  outputs  of 
interest,  [a,  q,  9,  n^,  y,  v],  while  Cp  relates  to  the  outputs  fed  back,  CL  and  q.  Input  weights,  W,,  and  W*, 
represent  sensor  noise  intensities  and  wind  gust  noise  intensities  respectively,  and  are  set  to  anticipated 
levels  and  altered  in  the  design  process  to  obtain  the  desired  noise  rejection  characteristics.  The  F  is  the 
second  column  of  A,  the  OC  column,  thereby  affecting  the  states  as  an  a  perturbation.  The  output  weight 

Wh,  is  a  diagonal  matrix  with  weighting  factors  on  each  of  the  outputs  of  H.  The  p  is  a  penalty  on  control 
power.  The  H2  problem  can  be  expressed  in  equation  form  as 


3-9 


From  this  set-up,  the  designer  sets  the  weights  based  on  expected  intensities  and  relative  importance  of  the 
variables,  then  tunes  the  weights  to  obtain  the  controller  that  exhibits  the  desired  objectives. 

One  of  the  desired  objectives  is  to  design  a  controller  so  the  closed  loop  system  will  track  a  given 
commanded  input.  Since  optimization  does  not  target  tracking  as  an  objective,  a  pseudo-integrator  was 
originally  placed  prior  to  the  actuator.  This  pseudo-integrator  would  be  joined  with  the  controller  for 
implementation  to  help  improve  tracking.  It  was  found  that  the  pseudo-integrator  acted  like  a  dynamic 
weight  on  u,  which  only  increased  the  order  of  the  controller,  complicated  the  numerics,  and  did  not 
improve  the  cracking.  The  H2  synthesis  created  a  controller  state  that  inverted  and  canceled  the  pseudo¬ 
integrator  state.  When  the  H2  controller  was  designed  without  the  pseudo-integrator,  the  steady  state 
response  could  be  scaled  with  an  input  gain  to  obtain  the  desired  outcome.  Therefore,  the  H2  controller  was 
designed  without  the  pseudo-integrator,  as  shown  in  Figure  3.5. 

The  process  of  designing  an  H2  controller  was  fairly  straight  forward  when  the  only  objective  was  to 
reject  noise.  Even  when  the  two  types  of  noise  introduced  into  the  system  resulted  in  competing  objectives, 
the  designer  could  quickly  determine  acceptable  levels  of  each  disturbance  and  trade  off  weightings  to  find 
an  acceptable  controller.  However,  designing  an  H2  controller,  based  on  this  structure,  that  met  the 
handling  qualities  criteria  was  extremely  difficult.  It  should  be  pointed  out  that  this  method,  as  well  as  the 
other  state-space  optimization  methods,  are  far  more  powerful  with  full  state  feedback.  The  process  of 
designing  an  H2  controller  was  just  as  iterative  as  the  classical  method,  but  more  complex.  For  the  classical 
method,  the  designer  chooses  a  solution  based  on  engineering  judgment,  evaluates  the  design  to  determine  if 
the  desired  objectives  were  met,  then  refines  the  design.  For  the  H2  controller,  the  designer  varies  the 
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problem  set-up  via  the  weightings,  solves  the  optimization  problem,  then  evaluates  the  design  to  determine 
if  the  desired  objectives  were  met.  The  designer  then  refines  the  weightings,  resolves  the  problem,  and 
hopes  the  ‘refined’  controller  is  an  improvement.  Unfortunately,  the  physical  interpretation  of  the  controller 
states  is  lost  in  the  optimization  process.  So,  although  choosing  which  weighting  to  change  to  increase  the 
noise  rejection  properties  on  a  given  output  variable  is  not  difficult,  knowing  which  weighting  would  effect 
the  closed  loop  damping  or  natural  frequency  was  not  intuitive.  There  was  no  predictable,  direct  relation 
between  output  weightings  and  improving  predicted  handling  qualities  ratings.  The  design  chosen  for 
implementation  was  found  after  numerous,  almost  random  iterations.  The  difficulty  in  achieving  the  desired 
handling  qualities  objective  by  this  design  method,  with  the  model  structure  set-up  used  in  this  research, 
raises  serious  questions  concerning  the  utility  of  this  design  method,  as  applied,  as  the  sole  source  for 
designing  a  manned  aircraft  flight  control  system.  Since  noise  rejection  typically  is  not  the  primary 
objective  of  a  flight  control  system,  H2  optimization  methods  are  well  suited  for  mixed  design  methods 
which  combine  other  objectives  with  noise  rejection. 

Despite  the  difficulties  in  obtaining  level  1  handling  qualities  using  H2  design  techniques  with  the 
set-up  applied  in  this  thesis,  a  controller  design  was  chosen  to  flight  test.  The  weightings  used  are  shown  in 
Table  3.3,  with  the  penalty  on  control  power  set  at  p=0.40. 


Table  3.3.  H2  Design  Weightings 


The  resulting  controller,  in  state-space  form  was 


-  53.52 

57.02 

0.00 

0.00 

2.86' 

57.24 

-  62.35 

-0.08 

0.00 

0.04 

31.54 

-31.40 

-0.05 

-0.56 

0.00 

-57.97 

61.66 

0.08 

0.00 

-0.04 

-  76.96 

-  210.86 

5.01 

-  43.33 

-  29.22 

(3.10) 
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-  56.24 

53.13 

61.50 

-56.24 

30.96 

-31.54 

-60.81 

57.97 

0.00 

0.00 

C2=[-3.85  -10.54  0.25  -2.17  -0.46] 

D2  =  [0.00  0.00] 


(3.11) 


(3.12) 

(3.13) 


and  in  transfer  function  form,  the  controller  can  be  expressed  in  short-hand  notation  as 

8,,  _  -292.46(20.00)(256)[0.28,  0.22] 
a  “(116.71)[0.74,  1952](- 0.33)(0.04) 

8,,  _  255.02(20.00)(-4.75)[0.30,  0.25] 
q  (116.71)[0.74,  19f  2](- 0.33)(0.04) 

where  (x)  represents  (s+x),  and  [^,  ct]  represents  (s^+2^as-KT^). 

The  time  histories  resulting  from  the  SIMULINK  simulation  are  shown  in  Appendix  C.  The 
handling  qualities  analysis  results  for  the  Bandwidth  criteria  are  shown  in  Figure  3.6.  and  the  results  for  the 
RSmith  criteria  are  shown  in  Table  3.4. 
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Parameter 

Value 

Level 

S  (dB/octave) 

-5.54 

1 

4.67 

-143.11 

2 

0.46 

1 

Table  3.4.  RSmith  Handling  Qualities  Predicted  Level  for  the  H2 


The  stability  vector  margins  for  the  H2  design  system  are  listed  in  Table  3.5. 


Margin  Type 

Vector  Gain  Margin 
(dB) 

Vector  Phase  Margin 
(degrees) 

Input  Margin 

[7.41  -5.32] 

±33.35 

Output  Margin 
{CL  loop  ‘broken’) 

IBBI 

±77.17 

Output  Margin 
(q  loop  ‘broken’) 

[7.92  -8.05] 

±35.16 

MIMO  Margins 
(sensitivity) 

[6.50  -3.68] 

±30.55 

MIMO  Margins 
(complementary  sensitivity) 

[3.27  -5.32] 

Table  3.5.  H2  Design  Vector  Stability  Margins 


3.4.3.  H„  Design  Set-up  and  Resulting  Design 

The  H„  design  problem  was  set-up  as  a  one  degree  of  freedom,  angle  of  attack  command  controller 
in  the  feed-forward  loop  with  dynamic  weightings  on  the  system  sensitivity  and  complementary  sensitivity 
and  a  weight  on  the  input  to  penalize  the  control  power  of  the  system.  Figure  3.7  depicts  the  block  diagram 
of  the  KL  problem  set-up. 


U 


63 


Figure  3.7.  BL  Design  Model  Block  Diagram 
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As  with  the  Hj  set-up,  the  actuator  is  separate  from  the  aircraft  plant,  and  Cp  contains  the  outputs  fed  back, 
OC  and  q.  The  input  weight,  Wj,  represents  a  weight  on  the  bounded  energy  exogenous  input,  d.  The 
dynamic  weights  on  S  and  T  are  simply  the  inverse  of  the  desired  loop  shapes  of  the  system’s  sensitivity  and 
complementary  sensitivity  functions.  The  p  is  a  penalty  on  control  power.  The  H„  problem  can  be 
expressed  in  equation  form  as 


■'act 

is 


/L 


BsCp 


Bp^act 

A 

^act 

0 

0 


Bu 


0 

0 

^s 

0 


0 

0 

0 

At'p 


. \ 

0 

0 

^act 

-t- 

0 

d  + 

Bact 

Bs^, 

0 

.Ur  . 

0 

0  _ 

9 

Deu 

^1 

/3. 

= 

'DsCp 

DpCp 

0 

0 

0  0 

0  0 

0  ■ 

Cr 

0 

^act 

.^T  . 

-t- 

'DsV/a 

0 

0 

d  + 

o' 

0 

.P. 

(3.16) 


(3.17) 


(3.18) 


The  process  of  designing  an  controller  is  similar  to  the  design  process  for  the  H2  controller.  The 
weightings  are  set  to  initial  values,  the  controller  is  solved  for,  then  the  resulting  controller  is  evaluated 
against  some  criteria.  Although  the  process  is  similar,  the  choice  of  weights  is  more  intuitive.  The  weights 
are  based  on  the  desired  loop  shape  of  the  system,  then  tweaked  to  obtain  the  desired  objectives.  Constant 
weights,  as  well  as  first  and  second  order  weights  on  sensitivity,  and  constant  and  first  order  weights  on 
complementary  sensitivity  were  extensively  investigated,  including  combinations  of  each  type  of  weight. 
Poles  in  the  sensitivity  weights  ranged  from  0.0001  -  10.0  radians,  and  the  zeros  investigated  ranged  from 
1.0  -  10,000  radians.  The  complementary  sensitivity  weights  considered  contained  poles  ranging  from  1.0  - 
10,000  radians,  and  zeros  ranging  from  0.0001  - 1.0  radians.  Based  on  this  problem  set-up,  with  the 
aerodynamic  model  used,  and  types  and  ranges  of  weights  considered,  the  best  resulting  design  for  handling 
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qualities,  time  response,  and  stability  margins  was  determined.  The  weights  chosen  for  sensitivity,  S,  and 
complementary  sensitivity,  T,  in  transfer  function  form  were 


^  _5(j+100)  (j+0.001) 

“  (j  +  0.001)  ~  (^+100) 


(5+ 0.001) 


The  resulting  controller,  in  state-space  form  was 


■-7.80e-l 

3.87e-l 

0 

6.05e-l 

-1.24e+0 

-8.04e-2 

Ule+O 

-7.20e-l 

-4.8k-2 

-6.09e-l 

2.40e-l 

8.04e-2 

-5.75e  +  2 

-2.79e  +  3 

6i3e-2 

0 

0 

0 

0 

0 

0 

-3.94e-3 

9.96e  - 1 

0 

1.78e-l 

1.78e-r 

llle-l 

1.77e-l 

-5.22e-l 

-5.24e-l 

2.73e-l 

2.74e-l 

0 

0 

4.45e-l 

0 

0 

4.49e-l 

1.76e-3 

1.77e-3 

0 

2.86e-f0 

0 

0 

157e-6 

0 

4.09e-2 

0 

0 

\57e-6 

-5.62e-l 

0 

0 

0 

—  4.62e  —  6 

0 

-4.09e-2 

0 

0 

2.42e-6 

-9.9Qe+2 

-6.20e+l 

-2.49-H3 

6.82e-^0 

4.78e-hl 

0 

0 

-LOOe-3 

0 

0 

0 

0 

0 

-l.OOe-3 

0 

0 

0 

0 

0 

-l.OOeH-2 

(3.19) 


(3.20) 


=[-6.40e-l-l  -3.11e-l-2  7.26e-l-l  -l.lOe  +  2  ^.67e-t-0  -2J7e  +  2  7.59e-l  5.32e  +  0]  (3.21) 


i)»=[0  0] 


(3.22) 


and  in  transfer  function  form,  the  controller  can  be  expressed  in  short-hand  notation  as 


8,*  _  -258.72(20.00)(153)[0.29,  0.24](0.19) 

(3.23) 

a  [0.75.  41.63](3.03)(-1.18)(0.064)(0.001) 

5,1,  -134.49(20.00)(133)[0.29,  0.24](-2.95e-7) 

(3.24) 

q  [0.75,  41.63](3.03)(-1.18)(0.064)(0.001) 
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The  time  histories  resulting  from  the  SIMULINK  simulation  are  shown  in  Appendix  C.  The 
handling  qualities  analysis  results  for  the  Bandwidth  criteria  are  shown  in  Figure  3.8.  and  the  results  for  the 
RSmith  criteria  are  shown  in  Table  3.6. 


Parameter 

Value 

Level 

S 

(dB/octave) 

-6.54 

1 

4.43 

HE3II 

(de.s) 

-155.11 

2 

^-peak 

(sec) 

0.58 

1 

Table  3.6.  RSmith  Handling  Qualities  Predicted  Level  for  the  H„  Design 


Figure  3.8.  Category  C  Bandwidth  Requirements:  KL  Design 


The  stability  vector  margins  for  the  HL,  design  system  are  listed  in  Table  3.7. 


Margin  Type 

Vector  Gain  Margin 
(dB) 

Vector  Phase  Margin 
(degrees) 

Input  Margin 

[  17.05  -7.60  ] 

±50.91 

Output  Margin 
(a  loop  ‘broken’) 

[  15.57  -16.78  ] 

±50.62 

Output  Margin 
(q  loop  ‘broken’) 

[  12.15  -oo  ] 

±68.78 

MIMO  Margins 
(sensitivity) 

[5.74  -3.42] 

±27.97 

MIMO  Margins 
(complementary  sensitivity) 

[3.26  -5.27] 

Table  3.7.  H„  Design  Vector  Stability  Margins 
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3.4.4.  Mixed  H2/H00  Design  Set-up  and  Resulting  Design 

The  purpose  of  the  mixed  design  technique  is  to  design  a  controller  that  optimally  trades  off  the 
objectives  of  both  the  H2  design  and  the  H«  design.  Upon  inspection  of  the  time  histories  of  the  H2  design 
and  the  fL,  design  discussed  in  previous  sections,  due  to  the  handling  qualities  objective,  there  is  little 
difference  between  the  noise  rejection  properties  of  the  two  designs.  Therefore,  determining  a  tradeoff 
between  these  two  designs,  as  proposed,  would  not  demonstrate  the  strengths  of  the  mixed  method.  The  H2 
portion  of  the  mixed  set-up  was  therefore  redone  to  pose  the  H2  problem  as  a  pure  noise  rejection  problem. 
The  nominal  H2  controller  used  in  the  mixed  problem  was  not  constrained  by  being  implementable  as  an 
independent  design,  that  is,  it  did  not  have  to  track  commands,  have  stability  margins,  nor  meet  any 
handling  qualities  criteria.  Although  applying  this  same  logic  to  the  mixed  design’s  IL  sub-problem  might 
result  in  a  mixed  H2/  H„  controller  with  the  best  tradeoff  between  stability  margins  and  noise  rejection, 
there  would  be  no  guarantee  that  any  of  the  resulting  controllers  would  meet  the  handling  qualities  criteria. 
The  Ho,  sub-problem  used  was  required  to  meet  both  the  stability  margin  requirements  and  the  handling 
qualities  objective,  but  no  restriction  on  noise  rejection  properties  or  other  implementation  issues  were 
made. 


Initially  the  BL  design  set-up  described  in  the  previous  section  was  used  for  the  mixed  H2/  IL 
synthesis.  Unfortunately,  the  complexity  of  a  single  H.  constraint  with  dynamic  weightings  on  both 
sensitivity  and  complementary  sensitivity  prevented  a  solution  to  the  mixed  control  problem  as  posed.  A 
new  KL  sub-problem  set-up,  only  weighting  the  system  complementary  sensitivity  function,  was  derived. 
The  resulting  BL  controller’s  handling  qualities  ratings  are  shown  in  Figure  3.10,  and  Table  3.9.  The  BL 
controller  expressed  in  shorthand  notation  is 

_  -2.07e-f-7(l000.00)(20.00)(l.53)[0.29,  0.24] 
oc  “  (2.11e-h5)[0.97,  762.38](2.92)(- 1.07)(0.060) 

_  -2.08e-t-7(l000.00)(20.00)(li:3)[0.29,  0.24] 
q  ~  (2.11e-h5)[0.97,  762.38](2.92)(- 1.07)(0.060) 
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This  design  was  not  implementable  due  to  the  second  order  high  frequency  poles  which  could  not  be 
eliminated  without  resulting  in  an  improper  controller  or  effecting  the  design  characteristics.  The 
implementation  issues  of  the  various  controllers  are  addressed  in  the  next  chapter. 

The  mixed  design  H2  sub-problem  used  the  same  structure  as  the  set-up  used  for  the  previous  H2 
design.  The  weights  chosen  for  the  mixed  design  H2  sub-problem  resulted  in  a  design  that  rejected  the  wind 
and  sensor  noise;  specifically  the  focus  was  on  minimizing  the  angle  of  attack  perturbations,  or  wind  gust 
noise.  The  weights  used  are  shown  in  Table  3.8. 


Value 

Value 

1000 

W„ 

200.00 

N„ 

1000 

Wn 

100.00 

_ Na _ 

1000 

Wf, 

0.01 

II 

1 - \ 

1 _ 1 

Wv 

0.01 

Wnz 

1.00 

Wv 

1.00 

1 

P 

1.0 

Table  3.8.  Design  Weighting  for  Mixed  H?  Sub-problem 


This  design  would  have  been  considered  a  Cooper-Harper  handling  qualities  rating  10  (uncontrollable). 
The  time  histories  for  this  design  are  shown  in  Appendix  C. 

The  mixed  H2/  HL  synthesis  was  performed  with  these  two  sub-problems.  The  initial  controller  used 
in  the  numerical  mixed  solution  is  the  H2  controller  described  in  the  previous  paragraph.  The  resulting  QL 
versus  y  curve  is  presented  in  Figure  3.9.  Each  of  the  points  on  the  curve  represent  a  controller  determined 
to  be  on  the  optimal  2-norm  versus  00-norm  curve.  Each  of  these  controllers  was  then  evaluated 


Figure  3.9.  Mixed  H2/BL  OC  versus  y  curve 
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for  handling  qualities,  stability  margins,  and  time  history  response.  There  did  not  seem  to  be  any 
predictable  pattern  to  the  handling  qualities  rating  resulting  from  the  closed  loop  systems  as  successive 
controllers  on  the  curve  were  analyzed.  A  few  had  acceptable  handling  qualities  predictions,  so  the  ‘best’ 
of  those,  considering  all  of  the  evaluation  criteria,  was  chosen  for  flight  test. 

The  resulting  mixed  controller,  in  state-space  form,  was 


r-5.00e  +  3 

-6.14e-l 

0 

0 

0 

0 

0‘ 

6.13e-l 

-1.50e  +  3  ■ 

-4.81e-2 
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0 

4.8  le- 2 

-8.35e-t-l 

251  e+Q 
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0 

0  - 

-251  e  +  Q 

-2.07e+l 

3A9e  +  \ 
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0 

(3.27) 
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-3.49e-t-l 

-2A2e  +  \ 

5.18e-)-0 
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0 

0 

0 

-5.18e-h0 

L62e  +  0 

-2.78e-2 

-  0 

0 

0 

0 

0 

2.78e-2 

-8i2e-2_ 

9f5e-2 

-3.71e-t-0' 

-3.01e-2 

1.05e-(-0 

-lOle-hl 

1.8  6e  +  2 

3.83e-t-l 

1.06e  +  2 

(3.28) 

6.46e  -1-  0 

1.66e  +  2 

-1.73e-Hl 

-2.02e  +  l 

_-5.83e-H0 

4.49e-t-l_ 

[-  1.50e  +  3 

-5.00e  +  3 

2.60e  +  0 

-1.48eH-0 

-3.68e-l 

3.54e-l 

-1.61e-2] 

(3.29) 

= 

[o  o] 

(3.30) 

and  in  transfer  function  form,  the  controller  can  be  expressed  in  short-hand  notation  as 

5,4  -139.0l[0.54,  2815.28][0.84,  40.96][0.48,  0.26] 

"V”  (4999.94)(1500.29)(83.44)[034,  41.65](-1.29)(0.085) 


5,4  581.82(-27623.10)[-0.093,  97.35](l9.83)[0.23,  0.26] 

~  ~  (4999.94)(1500.29)(83.44)[054,  41.65](- 1.29)(0.085) 


The  time  histories  resulting  from  the  SIMULINK  simulation  are  shown  in  Appendix  C.  The 
handling  qualities  analysis  results  for  the  Bandwidth  criteria  are  shown  in  Figure  3.10.  and  the  results  for 


the  RSmith  criteria  are  shown  in  Table  3.9. 


Sub-problem 

Mixed  H2/H..., 

Design 

Parameter 

Value  (Level) 

Value  (Level) 

S  (dB/octave) 

-5.79 

(1) 

-5.90  (1) 

4.61 

4.58 

-132.22 

(1) 

-146.58  (2) 

tq-pcak  (sec) 

0.041 

(1) 

0.48  (1) 

Table  3.  9.  RSmith  Handling  Qualities  Predicted  Level  for  the  Mixed  Design 
The  stability  vector  margins  for  the  mixed  design  system  are  listed  in  Table  3.10. 


HL  Sub-problem 

Mixed  H2/H„  Design 

Margin  Type 

VGM 

(dB) 

VPM 

(degrees) 

VGM 

(dB) 

VPM 

(degrees) 

Input  Margin 

[37.79  -11.78] 

±59.15 

[  10.77  -12.09] 

±44.13 

Output  Margin 
(a  loop  ‘broken’) 

■B 

±61.14 

IBB 

±75.96 

Output  Margin 
(q  loop  ‘broken’) 

msm 

±68.64 

[  11.04  -18.70] 

±52.45 

MIMO  Margins 
(sensitivity) 

[  7.25  -3.89] 

±42.05 

[  7.51  -3.97] 

±40.28 

MIMO  Margins 
(complementary  sensitivity) 

[  4.70  -10.89] 

[  4.55  -10.13] 

Table  3.10.  Mixed  Design  Vector  Stability  Margins 
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3.5.  Summary 

This  chapter  has  described  the  steps  taken  to  develop,  design,  and  evaluate  four  longitudinal  flight 
control  systems,  starting  with  stability  derivatives  describing  the  basic  longitudinal  modes  of  motion  of  a 
fictitious  unstable  fighter  type  aircraft.  The  methods  used  to  design  these  flight  control  systems  included 
classical  methods,  as  well  as  modem  state-space  methods,  specifically  H2,  HL,  and  mixed  H2/KL  control 
synthesis.  The  performance  of  each  design  was  not  intended  to  be  compared  to  the  other  designs.  The 
designs  offered  here  may  not  be  the  ‘best’  each  of  these  design  methods  can  create.  The  intent  was  merely 
to  demonstrate  the  ability  to  implement  these  design  methods  from  design  model  set-up  through  the  flight 
test  process.  The  practicality  of  these  design  methods  was  addressed  in  this  chapter,  and  will  be  further 
addressed  in  Chapter  4,  Implementation.  With  that  said,  the  designs  are  summarized  in  the  following  figure 


and  tables.  Time  history  plots  are  in  Appendix  C. 


Classical 

H2  Design 

Ho  Design 

Sub-problem 

Mixed  H2/H», 

Parameter 

Value  (Lvl) 

Value  (Lvl) 

Value  (Lvl) 

Value  (Lvl) 

Value  (Lvl) 

S  (dB/octave) 

-3.97  (1) 

-5.54  (1) 

-6.54  (1) 

-5.79  (1) 

-5.90  (1) 

5.05 

4.67 

4.43 

4.61 

4.58 

-150.84  (2) 

-143.11  (2) 

-155.11  (2) 

-132.22  (1) 

-146.58  (2) 

0.48  (1) 

0.46  (1) 

0.58  (1) 

0.041  (1) 

0.48  (1) 

Table  3.11.  RSmith  Handling  Qualities  Predicted  Level  Summary 


Figure  3.11.  Category  C  Bandwidth  Requirements:  All  Designs 
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Classical 

Ha  Design 

EL  Design 

Mixed 

Sub-problem 

Mixed  Ha/H^ 

Margin  Type 

[VGM]  (VPM) 
(dB)  (deg) 

[VGM]  (VPM) 
(dB)  (deg) 

[VGM]  (VPM) 
(dB)  (deg) 

[VGM]  (VPM) 
(dB)  (deg) 

[VGM]  (VPM) 
(dB)  (deg) 

[6  -oo]  (±61) 

[  7  -5  I  (±33) 

[17  -8]  (±51) 

[38  -12]  (±59) 

[11  -12]  (±44) 

Output  Margin 
(a  loop  ‘broken’) 

[12  -8]  (±43) 

[25  -cc]  (±77) 

[16  -17]  (±51) 

[28  -=]  (±61) 

[25  -~]  (±75) 

Output  Margin 

(q  loop  ‘broken’) 

[7  -9]  (±37) 

[8  -8]  (±35) 

[12  -~]  (±69) 

[13  -oo]  (±69) 

[11-18]  (±52) 

MIMO  Margins 

(sensitivity) 

[  3  -2]  (±18) 

[  6  -41  (±31) 

[  6  -31  (±28) 

[  7  -4]  (±42) 

[  8  -4]  (±40) 

MIMO  Margins 
(comp  sensitivity 

[  2  -31  (±18) 

[  3  -51  (±31) 

— 

[5-111  (±42) 

[  5  -101  (±40) 

Table  3.12.  Vector  Stability  Margin  Summary 
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IV.  Implementation 


4. 1.  Implementation  Model 

The  implementation  issues  raised  in  this  chapter  are  directed  at  specific  challenges  resulting  from  the 
limitations  of  the  Calspan  VSS  II  Learjet.  Some  of  the  concerns  are  general  in  nature  and  would  apply  to 
any  system  designed  with  these  methods,  while  other  concerns  are  specific  to  the  Calspan  VSS  11  Learjet 
system  and  the  flight  test  limitations  for  this  project.  This  chapter  addresses  the  factors  that  required  major 
attention  to  ensure  these  flight  control  systems  were  successfully  flight  tested.  The  first  problem  was  due  to 
the  VSS  II  simulation  computer  system’s  computing  speed  and  the  integration  routine  it  uses.  The  second 
issue  concerned  the  problems  that  result  when  implementing  an  unstable  controller.  The  last  major  problem 
was  due  to  the  actuator  model  used  in  the  design  process.  The  last  section  of  this  chapter  discusses  the 
model  verification  and  validation  testing  conducted  on  the  implementation  used  for  the  actual  flight  test. 

The  architecture  of  the  controller  implemented  needed  to  be  general  enough  to  accommodate  all  four 
designs.  This  requirement  would  allow  designs  to  be  quickly  and  efficiently  switched  in-flight,  saving 
valuable  flight  test  resources.  Four  separate  set-up  script  files  were  used  to  initialize  the  system  variables  at 
the  stroke  of  a  button. 

4. 1. 1.  VSS  II  System  Limitations 

The  implementation  requirements  for  VSS  11  provided  by  Calspan  are  presented  in  Chapter  II.  The 
requirements  that  effected  the  implementation  model  most  were:  the  integration  routine  (Euler)  and  step  size 
(O.Olsec)  used;  casting  the  controllers  in  pole/zero/gain  form  using  scalar  parameters  and  integrators,  like 
in  Figure  2.13;  and  ensuring  that  the  integrators  can  be  initialized  and  zero-ized  (reset)  in-flight  with  the 
system  running. 

The  integration  routine  used  by  the  VSS  11  computer  was  not  known  until  after  the  design  process 
was  complete.  It  is  recommended  that  the  design  phase  use  the  integration  routine  implemented  in  the 
aircraft.  Unfortunately,  that  was  not  known  during  the  design  phase.  Once  this  limitation  was  discovered, 
the  computer  simulations  were  redone  using  Euler  integration.  These  revised  time  histories  were 
numerically  unstable.  The  designs  that  had  problems  contained  poles  and/or  zeros  with  a  natural  frequency 
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greater  then  100  radians/second.  Dynamics  at  those  frequencies  are  beyond  the  range  which  affect  aircraft 
dynamic  modes  of  motion.  Therefore,  to  eliminate  the  implementation  problem,  those  poles  and  zeros  were 
removed  from  the  controller’s  transfer  functions.  The  system  norms,  stability  margins,  handling  qualities 
rating,  and  time  histories  were  reevaluated  to  ensure  this  change  did  not  significantly  effect  the  properties  of 
the  controller.  For  the  designs  implemented,  the  evaluation  criteria  were  unaffected  by  the  reduced  order 
controllers.  The  reduced  transfer  functions,  in  shorthand  notation,  chosen  for  implementation  are  shown  in 
Table  4. 1 .  Although  the  H2  controller  maintained  a  high  frequency  pole 


Alpha  Controller,  ^ 
a 

q  Controller,  ^ 

9 

Classical  Design 

34 

2.3 

(10) 

H2  Design 

-  292.46{20.00)(236)[0.28,  0.22] 

255.02(20.00)(-4.75)[0.30,  0.25] 

( 1 1 6.7 1)[0.74,  1 9.52|(-  0.33)(0.04) 

( 1 1 6.7 1)[0.74,  1 9.52](-  0.33)(0.04) 

H,„  Design 

-  258.72(20.00)(li3)[0.29,  0.24](0.I9) 

-  1 34.49(20.00)( I53)[0.29,  0.24](-  2.95e  -  7) 

[0.75,  41.63](3.03){-  1.I8)(0.064)(0.001) 

[0.75,  41.63](3.03)(-  1.18)(0.064)(0.00l) 

Mixed  H2/H„  Design 

-146.88f0.84.  40.96][0.48,  0.26] 

-2.14[0.84,  40.96][0.48,  0.26] 

(83.44)[0.54,  4r.65](-  I.29)(0.085) 

(83.44)[0i4,  41.65](-  1.29)(0.085) 

Table  4.1.  Reduced  Order  Controllers 

at  1 16.7  rad/sec,  which  could  have  been  discarded,  the  computer  SIMULINK  simulation  time  histories 
remained  stable  and  appeared  to  meet  the  Calspan  requirements. 

The  mixed  H„  sub-problem  design  was  eliminated  from  implementation  since  it  could  not  be  reduced 
to  a  lower  order.  It  contained  one  zero  and  three  poles  with  natural  frequencies  well  above  100  rad/sec. 
Removing  them  would  result  in  an  improper  system  or  significant  changes  to  the  properties  of  the 
controller.  This  limitation  did  not  prevent  this  controller  design  from  use  in  the  mixed  control  synthesis. 
However,  it  does  point  out  one  concern  of  these  methods,  as  applied.  Controllers  can  be  synthesized  which 
are  not  implementable  due  to  high  frequency  poles  or  zeros.  Consideration  should  be  given  to  adding  a 
constraint  to  the  problem  which  limits  the  frequency  of  the  resulting  controller. 

The  next  task,  in  developing  the  implementation  model,  was  to  cast  the  controller  model  in 
pole/zero/gain  form,  like  in  Figure  2.13.  This  single  controller  architecture  had  to  accommodate  all  four 
controller  designs,  even  though  each  controller  is  of  different  order,  and  the  classical  controller  is  in  the 
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feedback,  rather  that  the  feed-forward  path.  There  are  numerous  ways  to  construct  a  system  that  meets  the 
Calspan  requirements,  the  controller  model  initially  chosen  is  shown  in  Figure  4.1.  This  model  met  all  of 


Figure  4.1.  Initial  Implementation  Model 


the  requirements  Calspan  speciHed.  The  problems  that  occurred  did  not  become  apparent  until  the  ground 
testing  phase  of  the  project. 

4.1.2.  Final  H2,  and  Mixed  H/  Implementation  Model 

There  were  two  major  implementation  issues  discovered  during  the  ground  testing  phase.  The  first 
was  a  result  of  the  controllers  designed  by  the  modem  methods  being  unstable;  that  is,  they  contained  poles 
in  the  right  half  s-plane.  The  second  concern  resulted  from  simulating  the  first  order  design  model  actuator 
with  the  actual  aircraft  s  fourth  order  actuator.  The  latter  of  these  concerns  was  simply  resolved,  as 
described  at  the  end  of  this  section.  The  implementation  of  an  unstable  controller  was  more  difficult. 

When  the  controllers  were  implemented,  as  shown  in  Figure  4.1,  during  the  ground  testing,  after 
approximately  23  seconds  the  system  went  unstable  and  diverged.  Upon  examining  the  controller  integrator 
states,  the  individual  integrators  diverged  upon  activation  of  the  controller.  The  closed  loop  system 
remained  stable  since  the  controller  output  was  a  small  difference  of  large  numbers.  As  the  states  grew 
exponentially,  the  numerics  broke  down  and  the  system  eventually  went  unstable.  Therefore,  the  general 
architecture  needed  to  be  posed  in  a  way  that  kept  the  integration  states  from  growing  exponentially.  The 
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controllers  were  designed  using  state-space  methods,  so  a  state-space  implementation  architecture  was 
attempted. 

The  implementation  architecture  redesign  was  accomplished  by  using  a  similarity  transformation 
matrix  on  the  state  space  representation  of  the  reduced  order  controllers.  The  reduced  order  controllers  were 
transformed  into  observer  canonical  form.  The  SIMULINK  block  diagram  was  then  created  to  represent  a 
general  controller  system.  The  block  diagram  that  was  implemented  is  shown  in  Figure  4.2. 


Figure  4.2. 
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and  xl-x6  are  the  individual  states.  The  ‘type-C’  and  ‘type-H’  blocks  select  the  appropriate  feedback 
structure,  and  the  ‘feedback  sign’  block  changes  the  polarity  of  the  feedback  loop.  The  off  diagonal 
elements  of  the  Ac  matrix  were  left  in  variable  form  so  controllers  of  lower  order  could  be  accommodated. 
This  representation  also  meets  the  Calspan  requirements,  and  remains  numerically  stable.  The  Calspan 
engineers  added  ‘balance  and  hold’  blocks  and  replaced  the  1/s  blocks  with  ‘engage  integrators’.  These 
Calspan  created  blocks  are  simple  functions  that  allow  the  integrators  to  be  initialized  and  reset  by  the 
computer  operator  in-flight.  It  should  be  noted  that  this  is  just  one  way  to  solve  the  implementation 
problem  resulting  from  unstable  poles  in  the  controller.  It  may  not  be  the  ‘best’  solution.  Some  of  the 
states  still  gradually  increased  to  fairly  large  numbers,  but  did  not  grow  exponentially,  and  did  not  cause 
numerical  problems  for  the  duration  the  system  was  engaged  during  flight  test  (approximately  8-10 
minutes  per  pattern).  Although  the  unstable  controllers  were  successfully  implemented  for  this  research,  a 
flight  control  designer  should  approach  unstable  controllers  with  caution  since  the  instability  complicates 
implementation  issues.  It  should  be  emphasized  that  this  research  was  for  a  point  design  under  controlled 
conditions,  and  not  for  a  broad  aircraft  envelope  or  varying  conditions. 

The  last  major  implementation  issue  concerned  the  actuator  model  used  in  the  design.  After  the 
implementation  model  was  finalized,  the  ground  testing,  as  well  as  the  initial  verification  and  validation 
flight  test  revealed  unpredicted  oscillatory  motion  in  some  of  the  time  histories.  The  Calspan  VSS  II  system 
simulates  different  aircraft  dynamics  by  varying  computer  gains  to  modify  the  elevator  actuator  command. 
The  VSS  n  system’s  actuator  is  extremely  fast  (approximately  200  rad/sec)  and  can  be  modeled  as  a  fourth 
order  system.  Modifying  the  fourth  order  dynamics  down  to  first  order  dynamics  with  computer  gains  was 
not  considered  effective.  The  simulation  would  have  been  more  accurate  if  a  second  order  actuator  model 
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had  been  used  for  design.  The  VSS  ITs  actuator  was  considered  fast  enough  that  the  resulting  aircraft 


response  should  not  have  been  adversely  effected,  and  the  added  _  dynamics  would  accurately 

5*  +  20 

represent  the  aircraft  targeted  by  the  designs.  However,  when  the  _  dynamics  were  removed  the 

s  +  20 

unwanted  response  disappeared.  The  resulting  system  response  is  discussed  in  the  following  verification 
and  validation  testing  section. 

4.2.  Verification  and  Validation  Testing 

The  simulated  aircraft  dynamics  and  the  four  flight  control  designs  were  ground  tested  to  verify  the 
implementation  accuracy,  from  9  September  1996  through  20  September  1996.  During  this  time  the 
previous  mentioned  implementation  problems  were  addressed.  In  addition,  two  verification  flights  were 
flown  by  a  CALSPAN  safety  pilot  and  a  USAF  Test  Pilot  School  staff  pilot  prior  to  actual  flight  test  of  the 
designs.  Model  validation  testing  was  accomplished  on  the  four  flight  control  system  designs  while  flying 
to  and  from  Palmdale  Airport,  on  the  flight  test  sorties.  Only  one  design  was  examined  on  a  single  sortie. 
Heavy  weight  data  were  obtained  flying  to  Palmdale  Airport  and  light  weight  data  during  the  return  trip  to 
Edwards  AFB.  Programmed  Test  Inputs  (PTI)  provided  the  required  frequency  sweeps  and  step  inputs. 

The  frequency  sweeps  provided  data  for  Frequency  Response  Analysis  (FRA).  Data  were  recorded  using 
the  VSS  II  onboard  system  in  a  MATLAB™  format.  The  sampling  rate  was  50  Hz.  Recorded  data 
parameters  are  listed  in  Table  Dl,  in  Appendix  D.  All  model  validation  testing  occurred  between  4,000  and 
6,000  feet  MSL.  The  aircraft  was  flown  at  final  approach  speed  (125  to  135  BCAS)  with  gear  down  and  20- 
degree  flaps.  The  frequency  range  of  the  programmed  sweeps  was  from  .1  to  10  Hz.  The  amplitude  was 
0.25  degree.  Duration  was  25  seconds  which  provided  approximately  1250  data  points.  The  FRA  was 
performed  using  PC-FRA,  Version  1.01,  developed  by  High  Planes  Engineering.  Step  inputs  were  used  to 
compare  predicted  time  histories.  The  steps  were  1  inch  of  stick  deflection  in  magnitude  for  3  seconds 
duration.  Time  histories  and  bode  plots  were  then  compared  to  predicted  responses  simulated  in 
MATLAB™. 

Model  validation  analysis  was  performed  using  the  data  shown  in  Appendix  E.  Bode  plots  generated 
from  MATLAB™  and  flight  test  data  for  the  closed  loop  transfer  functions  alpha  output  to  stick 
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displacement  and  pitch  rate  output  to  stick  displacement  were  analyzed  for  each  design  method.  The  flight 
test  data  were  limited  in  frequency  range  from  0.1  radians  per  second  (r/s)  to  approximately  30  to  40  r/s, 
based  on  the  frequency  sweep  used.  The  higher  frequency  content  was  not  considered  significant  since  this 
was  beyond  the  frequency  range  that  affects  handling  qualities  and  most  of  the  system  outputs  in  the  flight 
control  systems  were  attenuated  above  30  r/s.  Lower  frequency  data  would  be  desired,  down  to  at  least  0.01 
r/s,  to  validate  the  low  frequency  response  of  the  closed-loop  system  and  to  compare  it  to  the  predicted 
model  since  there  were  some  dynamics  predicted  at  lower  frequencies.  Model  validation  testing  should  be 
conducted  with  frequency  sweeps  that  provide  a  lower  frequency  content. 

For  the  classical  flight  control  system  design,  the  magnitude  gain  over  the  frequency  range  tested 
matched  well  with  predicted  values.  The  H2  design  frequency  plots  showed  significant  differences  between 
the  predicted  gain  levels  and  the  flight  test  data.  The  H«.  design  frequency  plots  were  inconclusive  due  to 
scatter,  and  were  marginal  at  best.  The  mixed  design  frequency  plots  had  the  best  comparison  with  the 
predicted  model,  but  still  contained  some  additional  gain  above  10  r/s.  The  phase  content  of  the  closed-loop 
system,  for  all  the  designs,  consistently  showed  additional  phase  lag,  as  much  as  approximately  20  degrees 
from  predicted.  This  could  degrade  the  predicted  handling  qualities,  based  on  RSmith  criteria  [SG79],  as 
much  as  one  to  two  Cooper-Harper  ratings  [CH69]. 

Flight  test  time  histories  of  alpha  and  pitch  rate  (q)  resulting  from  a  step  programmed  test  input  (PTI) 
were  compared  to  MATLAB™  generated  model  predictions  (Figures  E9  through  E12).  The  predictions 
used  a  3-second  step  input,  while  the  duration  of  the  PTI  lasted  from  5  to  6  seconds.  Response  trend 
comparisons  could  still  be  made.  The  flight  test  pitch  rate  response  was  representative  of  the  predicted 
response  for  all  the  designs  except  the  H2  design,  which  deviated  from  predicted.  The  classical  and  the 
mixed  designs  exhibited  some  additional  pitch  rate  response  delay,  compared  to  the  predicted  model.  The 
alpha  flight  test  response  for  all  the  designs  was  similar  in  character  to  predicted  for  the  first  1  to  2  seconds, 
although  the  classical  design  displayed  more  initial  delay.  The  alpha  response  had  significant  deviations 
which  can  be  seen  in  Appendix  E,  Figures  E9  through  E12. 

The  close  matches  between  the  classical  predicted  curves  and  the  flight  test  data,  for  both  frequency 
response  and  time  histories,  tend  to  indicate  that  the  unaugmented  aerodynamics  simulated  by  the  VSS  II 
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Learjet  were  similar  to  the  model  used  in  design,  although  the  classical  system  exhibited  significant 
additional  lag.  This  in-flight  simulated  aerodynamic  model  was  used  for  all  the  configurations  flown.  The 
marginal  at  best  (mixed),  inconclusive  to  marginal  (IL.),  and  mismatched  (H2)  comparisons  between 
predicted  curves  and  flight  test  data  indicate  a  problem  in  the  implementation  of  these  designs  in  the  VSS  II 
Leaijet.  Although  there  is  a  mismatch  between  the  design  predicted  time  histories  and  the  model  validation 
flight  test  results,  for  the  purpose  of  demonstrating  the  design  process  through  flight  test  implementation  of 
these  methods,  the  difference  was  not  significant  enough  to  invalidate  the  results.  Analysis  of  the  VSS  II 
Learjet  system  capabilities  and  the  limitations  when  implementing  these  flight  control  designs  are  beyond 
the  scope  of  this  thesis,  but  should  be  considered  in  future  flight  tests  of  similar  systems.  It  should  also  be 
noted  that  these  factors  may  contribute  to  but  do  not  fully  explain  the  undesirable  aircraft  response  during 
the  flare  through  landing  roll-out,  exhibited  by  some  of  the  flight  control  systems  described  in  the  following 
results  sections. 

The  deviations  between  the  predicted  time  histories  and  the  flight  test  results  make  drawing 
conclusions  concerning  comparisons  between  design  methods  questionable.  This  was  not  an  objective  of 
this  research.  The  model  verification  and  validation  testing  results  were  considered  sufficient  to 
demonstrate  the  successful  implementation  of  these  flight  control  design  methods. 
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V.  Flight  Test  Results 


This  chapter  presents  the  results  of  a  limited  handling  qualities  evaluation  of  the  four  longitudinal 
flight  control  systems  designed  during  this  research.  Three  of  the  flight  control  systems  were  designed  using 
modem  state-space  optimization  techniques.  The  fourth  was  designed  using  classical  techniques.  The 
flight  test  phase  of  this  project  was  labeled  HAVE  INFINITY.  The  purpose  of  the  HAVE  INFINITY  flight 
test  was  to  evaluate  the  handling  qualities  of  these  four  longitudinal  flight  control  systems  during  the 
approach  and  landing  phase  of  flight.  Flight  tests  were  conducted  using  the  CALSPAN  Variable  Stability 
Simulator  II  (VSS  11)  aircraft  from  30  September  to  3  October  1996  at  Edwards  AFB  and  Palmdale  Airport, 
California.  The  HAVE  INFINITY  test  team  consisted  of  members  of  USAF  Test  Pilot  School  Class  96A. 
Team  members  flew  9  practice  sorties,  requiring  12  flight  hours,  and  6  test  sorties,  requiring  8  flight  hours. 
Testing  was  requested  by  Wright  Laboratory,  Wright-Patterson  AFB,  Ohio. 

5.1.  Test  Objective 

The  test  objective  was  to  evaluate  the  handling  qualities  of  the  four  HAVE  INFINITY  longitudinal 
flight  control  systems  during  the  approach  and  landing  phase  of  flight.  Each  flight  control  system  was 
evaluated  using  an  offset  landing  task  as  described  in  the  Test  Procedures  section  of  this  chapter.  Cooper- 
Harper  (CH)  and  PIO  ratings  were  assigned  by  the  project  test  pilots  after  each  offset  landing  task  . 

Analysis  of  CH  and  PIO  ratings,  as  well  as  pilot  comments,  allowed  the  determination  of  a  final  handling 
quality  level  for  each  flight  control  system  design.  One  of  three  levels  was  assigned  based  upon  the 
following  guidance. 

Level  I:  Satisfactory,  no  improvement  necessary,  CH  1-3 
Level  II:  Unsatisfactory  but  tolerable,  deficiencies  warrant  improvement,  CH  4-6 
Level  III:  Unacceptable,  deficiencies  require  improvement,  CH  7-9 
These  final  criteria  were  based  on  information  described  in  “The  Use  of  Pilot  Ratings  in  the 
Evaluation  of  Aircraft  Handling  Qualities”[CH69]  and  in  USAF  TPS  Flying  Qualities  Phase  Text  Vol.  IV, 
Flying  Qualities  Testing  [TPS95]. 
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5.2.  Test  Procedures 


The  HAVE  INFINITY  test  program  was  conducted  from  30  September  to  3  October  1996  at 
Edwards  AFB  and  Palmdale  Airport,  California.  It  consisted  of  six  flights  and  8  hours  of  flight  test.  The 
handling  qualities  of  four  different  longitudinal  flight  control  systems  were  evaluated  during  the  approach 
and  landing  phase  of  flight  using  an  offset  landing  task.  The  task  was  performed  by  four  evaluation  pilots 
with  a  broad  range  of  flying  experience.  Table  5.1  summarizes  each  pilot’s  previous  weapon  system 
experience. 


Evaluation 

Pilot 

Weapon  System  Experience 

1 

C-141B 

2 

B-IB 

3 

F-16 

4 

F-18 

Table  5.1.  Evaluation  Pilots’  Flying  Experience 


Each  of  the  four  evaluation  pilots  rated  the  flight  control  system  designs  using  the  CH  and  PIO  rating 
scales  (Appendix  F).  Prior  to  the  actual  flight  test,  each  pilot  practiced  the  offset  landing  task  in  a  variety  of 
different  aircraft.  This  allowed  the  pilots  to  become  familiar  with  the  task  over  a  broad  range  of  aircraft 
handling  qualities.  The  practice  aircraft  included  the  F-15,  F-16,  T-38,  and  C-23. 

5.2.1  Methods  and  Conditions 

The  test  team  aircrew  onboard  the  VSS  II  included  a  CALSPAN  safety  pilot,  a  HAVE  INFINITY 
evaluation  pilot,  a  test  conductor,  and  a  CALSPAN  system  operator  who  reconfigured  the  flight  control 
systems.  The  CALSPAN  safety  pilot  was  the  aircraft  commander.  The  test  conductor  duties  were 
performed  by  a  HAVE  INFINITY  test  team  engineer  or  a  second  evaluation  pilot.  Flight  tests  were  limited 
to  a  maximum  head-wind  of  20  knots,  a  maximum  tailwind  of  10  knots,  and  a  maximum  crosswind  of  10 
knots  for  flight  safety  and  data  quality  considerations.  The  ground  test  team  consisted  of  two  spotters 
stationed  beside  the  landing  runway  at  Palmdale  Airport.  These  spotters  maintained  VHF  radio  contact 
with  both  the  Palmdale  control  tower  and  the  test  team  aircrew. 
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The  flight  test  was  structured  into  test  blocks.  Each  test  block  represented  one  of  the  four  flight 
control  system  configurations  and  consisted  of  three  approaches.  The  three  approaches  had  to  be  flown 
sequentially  on  a  single  flight  in  order  to  complete  the  test  block.  Each  test  block  was  reserved  for  a  single 
pilot/flight  control  configuration  pairing.  Only  the  CALSPAN  safety  pilot  and  the  test  conductor  knew 
which  flight  control  system  was  actually  being  tested.  Evaluation  pilots  performed  blind  testing  throughout 
the  test  program. 

The  first  approach  of  each  test  block  was  a  straight-in  to  a  touch  and  go.  If  the  flight  control 
configuration  under  test  was  controllable,  as  defined  in  MIL-STD-1797A  [Mil90],  the  test  proceeded  to  the 
offset  landing  task.  The  offset  landings  were  accomplished  as  described  in  the  next  section.  The 
evaluation  pilot  provided  comments  and  CH/PIO  ratings  for  each  offset  landing  performed.  If  a  CH  rating 
of  8  or  better  was  obtained  on  the  first  offset  landing  of  a  test  block,  the  evaluation  pilot  repeated  the  offset 
landing  task  and  again  provide  comments  and  CH/PIO  ratings.  A  third  offset  landing  was  accomplished  if 
the  comments  and  ratings  from  the  first  two  were  inconsistent.  The  test  conductor  determined  if  the  ratings 
were  consistent  and  when  the  test  block  was  complete.  Data  from  any  additional  landings  were  also 
included  in  the  final  analysis  of  the  appropriate  control  system  design. 

5.2.2  Ojfset  Landing  Task: 

The  offset  landing  task  was  performed  at  Palmdale  Airport,  Runway  25.  The  VSS  II  aircraft  was 
configured  with  the  appropriate  flight  control  system  on  downwind  at  approximately  150  KCAS  and  the 
variable  stability  system  was  then  engaged.  The  aircraft  was  configured  with  gear  and  20-degree  flaps  and 
turned  to  final  to  intercept  the  3-degree  ILS  glideslope.  On  final,  the  aircraft  was  offset  300  feet  to  the  left 
of  centerline  and  slowed  to  final  approach  speed.  The  safety  pilot  called  altitude  above  the  ground  in  50- 
foot  increments  off  of  the  radar  altimeter.  At  200  feet  AGL,  the  safety  pilot  called  “maneuver”.  The 
evaluation  pilot  aggressively  maneuvered  the  aircraft  using  30  to  45  degrees  of  bank  to  capture  the  runway 
centerline.  A  simultaneous  longitudinal  correction  was  made  to  intercept  a  visual  glidepath  to  touchdown. 
The  pilot  attempted  to  land  in  the  middle  of  the  desired  box  described  in  the  next  paragraph  at  a  touchdown 
speed  10  knots  less  than  approach  speed  (+10/-5  knots). 
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5.2.2. 1.  Landing  Zone 

Specialized  runway  markings  used  to  determine  desired  and  adequate  task  performance  were  painted 
on  Runway  25  at  Palmdale  Airport.  Standard  18-inch  wide  white  paint  lines  were  used  for  all  markings. 

The  landing  zone  markings  are  depicted  in  Figure  5.1.  The  leading  edge  of  the  desired  box  corresponded  to 
the  ELS  glideslope  point  of  impact.  Standard  runway  distance  markers  were  used  as  a  backup  in  case  the 
lines  became  obscured  or  a  runway  change  was  necessary. 


600  Feet 


GRASS 


A  /  CONCRETE 


200  Feet 

400  Feet 

400  Feet 

Figure  5.1.  Landing  Zone  Markings  and  Dimensions 


5.2.2.2.  Landing  Task  Evaluation 

Following  the  touch  and  go,  the  safety  pilot  took  control  of  the  aircraft.  The  evaluation  pilot 
received  a  performance  rating  from  the  ground  spotters  over  the  VHF  radio  based  on  the  aircraft’s 
touchdown  position.  For  a  landing  to  be  considered  in  a  particular  zone,  both  main  landing  gear  were 
required  to  be  on  or  inside  the  respective  white  line.  The  evaluation  pilot  determined  whether  airspeed  and 
lateral  touchdown  criteria  were  met.  Based  on  this  information,  the  evaluation  pilot  assigned  a  performance 
rating  of  desired,  adequate,  or  inadequate.  The  test  conductor  then  stepped  the  evaluation  pilot  through  the 
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CH  and  PIO  rating  scales.  Evaluation  pilots  assigned  a  CH  and  PIO  rating  to  each  offset  landing  task 
accomplished.  Evaluation  pilots  also  gave  qualitative  comments  on  handling  qualities,  workload,  and 
problems  effecting  task  performance.  The  comment  card  used  is  shown  in  Appendix  F. 

An  8mm  camcorder  onboard  the  VSS  II  recorded  all  landings  and  associated  pilot  comments.  A 
ground  based  camera  also  recorded  Sortie  No.  2  from  final  approach  to  roll  out.  Post-flight  analysis  of 
cockpit  videotape  and  the  test  card  comments  was  used  to  construct  a  database  of  pilot  comments  for  each 
of  the  flight  control  system  designs.  A  summary  of  these  pilot  comments  can  be  found  in  Appendix  G. 

5.3.  Results  and  Analysis 

The  next  four  sub-sections  describe  the  flight  test  results  for  the  four  design  configurations  based  on 
the  flight  test  team’s  analysis  of  the  flight  test  data.  The  focus  of  the  HAVE  INFINITY  test  team  was  to 
determine  a  handling  qualities  rating  for  each  design.  The  test  team’s  perspective  and  evaluation  is 
included  even  though  it  was  not  limited  to  the  design  point  used  in  the  controllers’  design.  All  of  the  figures 
referenced  in  these  next  four  sub-sections  are  presented  in  Appendix  H,  and  are  as  they  were  presented  in 
the  flight  test  report.  The  overall  results  section,  however,  interprets  these  results  as  they  apply  to  this 
thesis. 

5.3.1.  Classical  Flight  Control  Configuration 

The  four  evaluation  pilots  performed  four  straight  in  approaches  and  ten  offset  landing  tasks  with  the 
classical  flight  control  design.  See  Figures  H.l.  and  H.2.  for  CH  and  PIO  ratings.  Evaluation  pilot  4  did 
not  provide  PIO  ratings  for  two  approaches.  All  the  pilots  stated  that  the  classical  flight  control  design  flew 
nicely  throughout  the  landing  task.  The  aircraft  was  predictable,  responsive  and  had  good  sensitivity.  One 
pilot  indicated  that  quick  abrupt  inputs  were  required  to  maintain  glide-path,  and  that  there  was  a  tendency 
to  over  flare  the  aircraft.  A  second  pilot  stated  that  he  got  a  slight  bobble  upon  the  initial  glide-slope 
capture  and  that  the  stick  was  a  bit  too  heavy.  The  third  pilot  stated  that  excessive  aft  stick  was  required  in 
the  flare.  A  fourth  pilot  stated  that  the  aircraft  flew  nicely  without  additional  qualifications.  Overall,  this 
configuration  was  considered  Level  II. 
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5.3.2.  H2  Flight  Control  Configuration 

The  four  evaluation  pilots  performed  four  straight  in  approaches  and  twelve  offset  landing  tasks  with 
the  H2  flight  control  design.  See  Figures  H.3.  and  H.4.  for  CH  and  PIO  ratings.  The  comments  indicated 
that  the  aircraft  felt  heavy  and  unresponsive.  One  pilot  called  the  aircraft  jittery.  The  aircraft  required  low 
gain  inputs  in  order  to  fly  acceptably.  Of  the  12  landings,  11  had  desired  touchdown  performance,  and  one 
was  inadequate.  All  pilots  indicated  that  the  aircraft  was  acceptably  flown  down  to  touchdown.  The 
aircraft  exhibited  an  uncommanded  nose  up  after  touch  down  requiring  extensive  forward  stick  input  to 
maintain  control.  On  one  occasion  the  aircraft  became  unexpectedly  airborne  after  touchdown.  The  aircraft 
performance  after  touchdown  drove  the  higher  CH  ratings.  Overall,  the  Ha  flight  control  configuration  was 
considered  to  have  Level  III  handling  qualities  due  to  the  unacceptable  motions  after  touchdown. 

5.3.3.  Flight  Control  Configuration 

The  four  evaluation  pilots  performed  five  straight  in  approaches  and  thirteen  offset  landing  tasks  with 
the  Hoo  flight  control  design.  See  Figures  H.5.  and  H.6.  for  CH  and  PIO  ratings.  The  comments  indicated 
that  the  aircraft  flew  predictably  until  the  pilot  increased  his  input  gains  when  the  aircraft  exhibited  a  jittery 
motion.  Of  the  13  landings,  6  had  desired  touchdown  performance,  3  were  adequate,  and  4  were 
inadequate.  Evaluation  pilot  4  did  not  provide  PIO  ratings  for  two  approaches.  One  pilot  noticed  a 
lightening  in  control  forces  with  reduction  in  airspeed  for  all  four  of  his  landings.  For  all  pilots  the  aircraft 
behaved  unpredictably  upon  flare  initiation.  The  aircraft  tended  to  float  and  all  pilots  noted  the  need  to 
apply  forward  stick  right  after  the  initiation  of  the  flare  to  arrest  the  pitch-up  motion.  After  the  main  gear 
touchdown,  the  nose  could  not  be  lowered  even  with  full  forward  stick  input.  On  at  least  one  occasion  the 
aircraft  unintentionally  ballooned  and  became  airborne  after  touchdown.  The  aircraft  performance  after  the 
initiation  of  the  flare  through  landing  roll  out  resulted  in  the  high  CH  ratings.  Overall,  the  EL  flight  control 
configuration  was  considered  to  have  Level  III  handling  qualities  due  to  the  unacceptable  motions  after  the 
flare  initiation  through  landing  roll-out. 
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5.3.4.  Mixed  H2/H„  Flight  Control  Configuration 

The  four  evaluation  pilots  performed  a  total  of  four  straight  in  approaches  and  ten  offset  landings 
with  the  mixed  flight  control  design.  See  Figures  H.7.  and  H.8.  for  CH  and  PIO  ratings.  The  comments 
indicated  that  the  aircraft  felt  sluggish  and  heavy  on  approach  with  a  somewhat  slow  initial  response.  Two 
of  the  pilots  indicated  that  the  aircraft  was  somewhat  unpredictable,  while  the  other  two  considered  the 
aircraft  predictable.  All  four  pilots  stated  that  the  aircraft  handling  qualities  prior  to  the  initiation  of  the 
flare  were  acceptable.  For  the  10  offset  landings,  6  resulted  in  desired  touchdown  performance,  2  resulted 
in  adequate,  and  2  were  inadequate.  The  flight  control  configuration  had  unacceptable  behavior  from  flare 
initiation  through  landing  roll-out.  All  pilots  commented  on  a  heave  and  float  tendencies,  and  on  at  least 
four  occasions,  the  aircraft  became  airborne  after  the  initial  touchdown.  On  every  landing,  all  pilots 
required  one-half  to  full  forward  stick  to  arrest  an  unanticipated  pitch-up  motion  after  touchdown.  Overall, 
the  mixed  flight  control  configuration  was  considered  to  have  Level  III  handling  qualities  due  to  the 
unacceptable  motions  after  flare  through  landing  roll. 

5.3.5.  Overall  Results: 

The  HAVE  INFINITY  flight  test  team  concluded  that  the  classical  longitudinal  flight  control  design 
exhibited  Level  11  handling  qualities  for  the  approach  and  landing  phase  of  flight.  Furthermore,  all  three 
state-space  designs  exhibited  undesirable  motion  at  some  point  between  flare  and  the  landing  roll,  resulting 
in  Level  IE  ratings.  However,  the  landing  roll-out  was  not  included  in  the  design  phase.  The  changes  in 
stability  derivatives  due  to  ground  effect  and  airspeed  reduction  during  landing  roll;  the  effects  of  a  sensed 
pitch  rate  change  due  to  the  gear  impacting  the  ground;  the  shift  in  the  center  of  rotation  once  the  gear  are  in 
contact  with  the  ground;  and  the  thrust  effects  on  pitching  moment  as  the  throttles  are  retarded  to  idle  are 
some  of  the  factors  which  were  not  included  in  the  design  point  used  to  define  the  controllers  synthesized  in 
this  thesis.  While  it  is  interesting  that  the  classical  controller  did  not  exhibit  the  undesirable  characteristics 
on  landing  roll  and  the  modem  controllers  did,  the  landing  roll  handling  qualities  were  not  the  target 
objectives  of  the  designs. 

The  flight  control  designs  were  solely  designed  for  the  approach  phase  of  flight  and  intended  to  be 

analyzed  from  approach  to  the  touchdown  point.  The  pilot  ratings  that  were  given  were  mostly  based  on  the 
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disturbing  handling  characteristics  after  touchdown.  Therefore  the  ratings  given  were  not  directly 
applicable  to  the  desired  evaluation  objective.  Pilot  comments  did  indicate  the  three  state-space  designs 
exhibited  acceptable  handling  qualities  on  final  approach  although  no  precise  task  was  defined  for  this 
portion  of  the  approach.  Additionally,  many  of  the  spot  landings  met  the  desired  objective  criteria.  While 
it  is  not  appropriate  to  alter  the  flight  test  handling  qualities  ratings  after  the  fact,  the  data  does  support  the 
conclusion  that  all  designs  were  ‘acceptable’  on  approach.  The  term  ‘acceptable’  must  be  used  loosely 
since  there  is  insufficient  data  (i.e.  no  pilot  ratings  given  just  for  the  approach  portion  of  the  flight)  to 
determine  a  handling  qualities  level  for  the  approach  flight  phase  separate  from  the  landing  roll.  It  can  be 
concluded  that  the  four  flight  controller  designs  were  successfully  implemented,  safely  flown  during 
approach  to  the  touchdown  point,  and  although  the  handing  qualities  level  rating  data  was  inconclusive  for 
the  objective  flight  condition,  the  pilots  did  not  find  the  handling  qualities  objectionable  during  the 
approach  phase  of  flight. 
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VI.  Conclusions  and  Recommendations 


6.1.  Summary 

The  main  objective  of  this  thesis  was  to  apply  modem  control  theory,  specifically  H2,  H«,  and  mixed 
H2/IL,  control  synthesis,  to  an  actual  aircraft  system,  starting  from  an  aircraft’s  stability  derivatives  and 
progressing  through  implementation  to  flight  test  of  the  systems.  That  objective  was  successfully  met. 

The  stability  derivatives  of  a  fictitious  aircraft  with  unstable  airframe  aerodynamics  was  the  starting 
point  for  this  research.  Those  stability  derivatives  were  used  to  create  an  aircraft  model  in  both  transfer 
function  and  state-space  form.  Four  longitudinal  flight  control  system  controllers  were  design  for  this 
fictitious  aircraft.  These  designs  applied  not  only  classical  design  methods,  but  also  optimal  flight  control 
design  methods,  specifically  H2,  HL,  and  mixed  H2/IL  control  synthesis.  The  TPS  Flight  Control  System 
project  was  the  basis  for  the  structure  used  to  design  the  feedback  control  systems.  These  designs  were 
analyzed  for  handling  qualities,  noise  rejection  properties,  and  stability  margins.  Controller  designs  were 
obtained  that  met  the  design  analysis  criteria. 

The  classical  controller  was  primarily  included  as  a  control  case  for  comparison  and  to  help  verify 
proper  implementation  of  the  un-augmented  aerodynamics  simulated  by  the  VSS  II  system.  This  base-line 
design  applied  the  design  techniques  commonly  used  to  design  flight  control  systems,  including  root  locus 
and  bode  plot  analysis,  to  synthesize  the  controller.  This  method  was  tedious,  but  the  physical 
interpretation  of  the  controller  states  made  the  process  easy  to  analyze  and  improving  the  predicted  handling 
qualities  ratings  was  intuitive.  There  was  no  way,  except  trial  and  error,  to  optimize  the  controller  for  noise 
rejection  or  stability  margins,  or  to  quantify  trade-offs  between  performance  objectives. 

The  H2  design  problem  was  set-up  as  a  one  degree  of  freedom,  angle  of  attack  conunand  controller  in 
the  feed-forward  loop.  The  H2  design  technique  was  very  capable  at  both  wind  gust  and  measurement  noise 
rejection,  but  difficult  to  use  when  handling  qualities  was  an  objective.  The  use  of  a  pseudo-integrator  to 
improve  tracking  was  not  effective.  Setting  up  the  system  as  a  model  matching  problem,  and  investigating 
the  closed  loop  system  pole  locations  as  a  function  of  the  design  weightings,  might  provide  insights  to  the 
missing  physical  interpretation  of  the  controller  states  and  make  designing  for  handling  qualities  more 
intuitive. 


6-1 


The  design  problem  was  set-up  as  a  one  degree  of  freedom,  angle  of  attack  command  controller 
in  the  feed-forward  loop.  This  method  was  less  cumbersome  then  the  H2  technique,  but  also  fails  to  provide 
physical  interpretation  to  the  controller  states.  Systematic  changes  in  weightings  does  not  systematically 
improve  the  handling  qualities  predictions. 

The  mixed  H2/H=o  control  synthesis  optimally  traded  off  properties  of  an  H2  controller  solely 
designed  for  noise  rejection  and  an  IL,  controller  designed  to  obtain  both  good  margins  and  reasonable 
handling  qualities  with  a  weighting  only  on  complementary  sensitivity.  Again,  the  resulting  mixed 
controller  lacked  physical  interpretations  of  the  controller  states  and  any  correlation  between  handling 
qualities  predictions  and  controller  selection. 

The  four  flight  control  designs  were  then  implemented  in  the  Calspan  Variable  Stability  Simulator  II 
(VSS  II)  Learjet.  This  aircraft  has  the  capability  to  simulate  the  fictitious  aircraft’s  aerodynamics  and 
implement  the  controller  designs.  The  modern  designs  were  difficult  to  implement  since  those  methods 
tended  to  create  controllers  with  high  frequency  poles,  which  were  not  realistically  meaningful.  Also,  the 
resulting  controllers  were  unstable,  and  although  unstable  controllers  can  and  were  implemented,  internal 
stability  can  become  a  problem.  Unstable  controllers  are  not  ideal  in  applied  systems.  The  implementation 
of  the  fictitious  aircraft  aerodynamics  appeared  correct,  based  on  the  classical  design  verification  testing 
results.  The  verification  and  validation  testing  of  the  closed  loop  systems  indicated  significant  variations 
between  the  predicted  flight  control  system  response  and  the  flight  test  results.  This  would  limit 
conclusions  comparing  the  design  methods  and  the  resulting  controllers,  since  the  systems  tested  did  not 
accurately  represent  the  flight  control  systems  designed.  However,  such  comparisons  were  not  a  thesis 
objective.  The  variations  in  predicted  versus  flight  test  results  were  not  significant  enough  to  invalidate  the 
implementation  of  these  designs.  The  modern  designs  were  intended  to  accommodate  uncertainty,  and  this 
was  effectively  demonstrated.  A  better  understanding  of  the  VSS  II  system  capabilities  and  limitations 
would  improve  future  flight  testing  results  utilizing  this  system. 

The  four  controllers  were  evaluated  during  a  limited  handling  qualities  flight  test,  named  HAVE 
INFINITY.  The  classical  design  was  rated  level  11,  while  all  three  modern  designs  were  rated  level  El.  The 
resulting  pilot  ratings  were  influenced  by  the  handling  characteristics  after  touchdown  through  landing  roll. 
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This  was  not  part  of  the  flight  condition  the  controllers  were  designed  to  handle.  Pilot  comments  regarding 
the  handling  qualities  prior  to  touch  down  indicate  ‘acceptable’  handling  qualities  for  all  designs  tested. 

The  flight  control  designs  were  safely  implemented  and  flown,  but  the  flight  test  was  inconclusive  in 
determining  the  handling  qualities  rating  for  the  controller  in  the  intended  flight  condition. 

This  study  encompassed  aerodynamic  modeling,  aircraft  equations  of  motion  and  dynamics,  linear 
control  theory,  optimal  flight  control  design  methods,  aircraft  stability  and  control  analysis,  handling 
qualities  analysis,  principles  of  flight  test,  and  flight  test  project  management.  In  the  process  of  putting  all 
of  these  disciplines  together,  issues  were  raised  and  areas  of  further  research  were  discovered.  This  thesis 
was  accomplished  in  three  major  phases:  Set-Up  and  Design;  Implementation;  and  Flight  Test  and  Analysis. 
The  remainder  of  this  chapter  will  discuss  recommendations  arising  from  each  of  these  phases. 

6.2.  Recommendations 

•  Examine  and  flight  test  these  design  methods  using  full  state  feedback.  The  controllers  designed 
in  this  thesis  were  limited  by  feeding  back  only  angle  of  attack  and  pitch  rate.  Full  state 
feedback  does  not  complicate  the  optimization  process,  but  would  improve  control  of  the  states. 

•  Set-up  the  system  as  a  model  matching  problem. .  First  choose  an  aircraft  model  that  exhibits 
solid  level  1  handling  qualities.  Then  apply  H2  synthesis  to  design  a  controller  that  optimizes 
noise  rejection  while  minimizing  the  difference  between  the  system  output  and  the  ‘level  1  ’ 
aircraft  model.  This  set-up  was  not  used  during  this  thesis,  but  might  prove  to  be  a  more 
efficient  approach  to  the  problem,  should  noise  rejection  be  a  primary  objective  of  the  flight 
control  design 

•  Investigate  how  the  closed  loop  system  poles  are  affected  by  variations  in  design  weights.  The 
physical  interpretation  of  controller  states  and  the  effects  of  weighting  changes  on  handling 
qualities  predictions  could  be  explored  by  a  systematic  investigation  of  the  effects  on  pole 
locations  due  to  variations  in  design  weights.  Creating  a  ‘root  locus’  of  the  closed  loop  poles  as 
a  function  of  design  weighting  might  provide  insights  into  these  areas. 
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•  Include  problem  constraints  in  the  mixed  design  which  limit  the  controller’s  poles  to  a  realistic 
frequency  range.  Some  acceptable  controllers  were  eliminated  since  they  were  not 
implementable  due  to  high  frequency  poles.  The  mixed  design  process  is  capable  of  multiple 
constraints  of  various  types. 

•  Investigate  constraints  which  would  define  the  handling  qualities  criteria  boundaries  for  use  with 
mixed  control  synthesis.  If  the  handling  qualities  criteria  could  be  part  of  the  constraints  on  the 
mixed  problem,  the  handling  qualities  issue  would  be  eliminated  form  the  process. 

•  Use  a  known  design  method  as  a  control  case,  like  the  classical  design,  to  verify  and  validate  the 
implementation  and  help  isolate  contributing  factors  of  the  experiment.  The  classical  controller 
also  aided  in  the  blind  evaluations  during  the  flight  test. 

•  Either  design  the  controllers  to  include  landing  roll-out  as  a  design  criteria,  or  limit  the  flight  test 
evaluation  to  the  flight  condition  used  in  the  design  process.  The  point  design  used  in  the  design 
phase  of  this  research  did  not  include  landing  roll  characteristics,  but  the  pilots  included  this 
phase  of  flight  in  the  ratings  given.  As  a  result  the  ratings  did  not  reflect  the  condition  the 
controllers  were  designed  to  accommodate.  The  only  applicable  flight  test  data  were  general 
pilot  comments  concerning  the  approach  phase,  insufficient  for  determining  a  handling  qualities 
level. 
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Appendix  A 
Aircraft  Model  Data 


The  following  was  provided  by  Calspan: 


AIRCRAFT  DATA 
Flight  Condition 


Altitude 

1000 

ft 

TAS 

=  214.2 

ft/sec 

Mach 

=  0.1926 

KCAS 

=: 

125 

knots 

Oo 

=  6.69 

degrees 

Aircraft  Geometry 

S 

= 

232 

ft^ 

b 

=  34.1 

ft 

cbar 

=  7.0 

Mass  Properties 

Weight 

= 

11,000 

lbs 

Ixx 

=28,000 

slug-ft 

lyy 

=22,000 

Izz 

= 

47,000 

slug-ft^ 

Longitudinal  Stability  Derivatives 

Cdo 

= 

0.0486 

Clo 

=  0.36 

CmO 

=  -  0.004 

Cdi 

= 

0.0792 

Clu 

=  0.0 

Cmu 

=  0.0 

Cdu 

=s 

0.0 

Clb 

=  5.04 

Cma 

=  0.2 

Ctxu 

-0.3255 

Ci_ad 

=  0.0 

Cmad 

=  0.0 

Code 

= 

0.0 

Ccq 

=  0.0 

Cmq 

=  -6.0 

Coig 

=z 

0.02 

Ctde 

=  0.2438 

Cmde 

=  -0.7314 

The  following  is  a  MATLAB''''^  script  file  which  uses  the  above  aircraft  data  and  calculates  the  state- 

space  representation  of  the  aircraft  dynamics. 

function  [A,B,C,D]=  cspanmod 
%  function  [A,B,C,D]  =cspanmod 
%  Returns  a  4th-order  wind  axis  state-space  model 
%  of  a  Learjet  N1 02VS  @  1 25  Ki AS  and  1 000  feet 
%  The  states  (x)  are: 


%  1  q 

deg/sec 

%  2  aipha 

deg 

%  3  V  (true) 

feeVsec 

%  4  gamma 

deg 

%  Inputs  (u)  are: 

%  1  de 

deg 

%  2  dz 

deg 

%  3  dx 

deg 

%  Outputs  (y)  are: 

%  1  q 

deg/sec 

%  2  alpha 

deg 

%  3  V  (true) 

feet/sec 

%  4  gamma 

deg 

%  5  q_dot 

deg/sec/sec 

%  6  alpha  dot 

deg/sec 

%  7V_dot 

feet/sec/sec 

%  9  Nz 

g's 

%  1 0  theta 

deg 

%  Constants 

9=32.2; 

%ft/sec 

r2d=1 80/pi: 

%deg/rad 

d2r=pi/180; 

%rad/deg 

%  Right  condition 

VJndicated=125; 

%  KIAS 

h=1000; 

%ft 

A-1 


%  Mass  properties 
lyy=22000: 
weight=1 1 ,000; 
m=weight/32.2; 
%  Geometry 
b=34.l: 

3=232; 

c=7; 

tailJo_cg=21 ; 

ar=b/c; 

e=0.84; 


%  (slug  feet''2) 
%  lbs 
%  slugs 

%ft 

%ft''2 

%ft 

%ft 


%  Atmosphere  model 

rho_0=0.0023769;  %  slug/ft''3 

rho=rho_0*(1  -(.001 981 2*h/288. 1  5))a((1  /(.001 9812*96.034))-1 ); 
slug/ft''3 

V=VJndicated/sqrt(rho/rho_0)*6080/3600;  %  feet  per  sec 

q_bar=1/2  *  rho 

%  Non-dimensional  derivatives  (per  radian) 
cm0=  -0.004; 
cma=  0.2; 
cmde=-0.7314; 
cmad=  0; 
cmq=  -6; 
cm=  0; 


clO=  0.36; 
cla=  5.04; 

clde=-cmde*c/taiLto_cg; 
cl=weight/(q_bar*S); 
alphaO=(cl-clO)/cla; 
cdO=  0.0486; 

cda=(2*cl0*cla+2*cla*cla*alpha0)/(pi*ar*e); 
cdi=  0.0792*r2d; 
cdlg=  0.02; 

cd=cdO+cda*alphaO+cdlg; 

%  Dimensional  derivatives 

Mq  =rho*S*V*c*c/(4'lyy)'(cmq+cmad); 
M_alpha=rho*S*V*V*c/(2*lyy)*cma; 

Mu  =0.0; 

Zq  =0.0; 

Zw  =rho*S*V/(2*m)*(-cla); 

Z_alpha=Zw*V; 

Zu  =-2*g/V; 

Xq  =0.0; 

Xw  =  rho*S*V/(2*m)*(-cda); 
X_alpha=Xw*V; 

Xu  =  rho*S*V/m*(-cd); 

M_de  =rho'S*V*V*c/(2*lyy)*cmde; 

Z_de  =rho*S*V*V/(2*m)*(-clde); 

X_de  =0; 


A=[ 

Mq  M  alpha 

Mu 

l-hZqA/  Z  alpha/V 

Zu/V*57.3 

0  X_alpha/57.3 

Xu 

0  -Z_alphaA/ 

-ZuA/*57.3 

B  =  [ 

M  de 

Z  de/V 

X  de 

Z_de/V  ]; 

C=[ 

eye(size(A)) 

A 

A(4,;)*(V*d2r/g) 

0101]; 

D  =  [ 

zeros(size(B)) 

B 

B(4,:)*(V*d2r/g) 

0 

0 


-g/57.3 

0 


% 
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Appendix  B 

Design  Evaluation  Script  Files 

The  following  is  a  MATLAB™  script  file  which  calculates  the  vector  stability  margins  given  a  state- 
space  representation  of  the  system  [Asys,Busys,Cysys,Dysys]  and  the  controller  [Ak,Bk,Ck,Dk]. 


% 

%  Margins  Tinput  vector  margins 

%  2-output  (alfa  loop  broken)  vector  margins 

%  3-output  (q  loop  broken)  vector  margins 

%  4-muiti  variabie  output  vector  margins 

[akg,bkg,ckg,dkg]=series(Asys,Busys,Cysys,Dysys,Ak,Bk,Ck,Dk): 

[akgs,bkgs,ckgs,dkgs]=feedbk(akg,-bkg,ckg,-dkg,1); 

iakgt,bkgt,ckgt,dkgt]=feedbk(akg,-bkg,ckg,-dkg,2); 

isvs1,ws1]=sigma(akgs,bkgs,ckgs,dkgs); 

[svt1  ,wt1  ]=sigma(akgt,bkgt,ckgt,dkgt); 
alphao1=1/max(svs1 ); 
ao1=1/max(svt1); 
vgmst  lo=(1/(1  -t-alphaol )); 
vgmst  hi=(1  /(1  -alphaol )); 

if  vgmst  hi  <  1e-10;  vgmst  hi=1e-10;end; 
vgmttlo=(t-aot): 

if  vgmttio  <  te-tO;  vgmttlo=te-tO;end; 
vgmtthi=(t-Haot); 

vgmtlo=min([vgmstlo  vgmst  hi  vgmttio  vgmtthi]); 
vgmthi=max([vgmstio  vgmsthi  vgmttio  vgmtthi]); 
vgmt=20’logt0([vgmthi  vgmtio]); 
vpmst  =2*asin(alphaot  /2)*t  80/pi; 
vpmtt  =2*asin(aot /2)*t  80/pi; 
vpmt=max([vpmst  vpmttj); 

fprintf('a:\margins.m','\nlnput  VGM  =  [  %4.2f  %4.2f  ]dB\n‘,vgmt); 
fprintf('a:\margins.m', 'Input  VPM  =  %4.2f  deg\n',vpm1); 

[agk,bgk,cgk,dgk]=series(Ak,Bk,Ck,Dk,Asys,Busys.Cysys,Dysys); 

acl2=(agk-t-bgk(:,2)*cgk(2,;)); 

bcl2=bgk(:,t); 

ccl2=cgk(t,:); 

dcl2=dgk(t,t); 

[agk2s,bgk2s,cgk2s,dgk2s]=feedbk(acl2,-bcl2,ccl2,-dcl2,t);  %feedbk  defined  for 
[agk2t,bgk2t,cgk2t,dgk2t]=feedbk(acl2,-bcl2,ccl2,-dcl2.2);  %  negitive  feedback 
isvs2,ws2]=sigma(agk2s,bgk2s,cgk2s,dgk2s); 
[svt2,wt2]=sigma(agk2t,bgk2t,cgk2t,dgk2t); 
alphao2=t/max{svs2); 
ao2=t/max(svt2); 
vgms2lo=(t/(t-Falphao2)); 
vgms2hi={t /(t -alphao2)); 
if  vgms2hi  <  te-tO;  vgmts2hi=te-t0;end; 
vgmt2lo=(t-ao2); 

if  vgmt2to  <  1e-t0;  vgmt2lo=te-t0;end; 
vgmt2hi=(t-hao2); 

vgm2lo=min([vgms2lo  vgms2hi  vgmt2lo  vgmt2hi]); 
vgm2hi=max([vgms2lo  vgms2hi  vgmt2lo  vgmt2hil); 
vgm2=20*logt0([vgm2hi  vgm2lo]); 
vpms2=2*asin(alphao2/2)*t  80/pi; 
vpmt2=2*asin(ao2/2)*t  80/pi; 
vpm2=max([vpms2  vpmt2]); 

acl3=(agk-)-bgk(:,t)*cgk(t,;)); 

bcl3=bgk(:,2); 

ccl3=cgk(2,;); 

dcl3=dgk(2,2); 

[agk3s,bgk3s,cgk3s,dgk3s]=feedbk(acl3,-bcl3,ccl3,-dcl3,1);  %feedbk  defined  for 
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[agk3t,bgk3t,cgk3t,dgk3t]=feedbk(acl3,-bcl3,ccl3,-dcl3,2):  %  negitive  feedback 
[svs3,ws3]=sigma(agk3s,bgk3s,cgk3s,dgk3s); 

[svt3,wt3]=sigma(agk3t,bgk31,cgk3t,dgk3t): 

alphao3=1/max(svs3): 

ao3=1/max(svt3); 

vgnns3Io=(1/{1+alphao3)); 

vgms3hi=(1/(1-alphao3)); 

if  vgms3hi  <  1e-10;  vgms3hi=1e-10;end: 
vgmt3lo=(1-ao3): 

if  vgmt3lo  <  1e-10;  vgmt3lo=1e-10;end; 
vgmt3hi=(1+ao3): 

vgm3lo=min([vgms3lo  vgms3hi  vgmt3lo  vgmt3hil); 
vgm3hi=max([vgms3lo  vgms3hi  vgmtSlo  vgmtShi]); 
vgm3=20*log10([vgm3hi  vgm3lo]): 
vpms3=2*asin(alphao3/2)*1 80/pi; 
vpmt3=2*asin(ao3/2)*1 80/pi; 
vpm3=max([vpms3  vpmt3]); 

fprintf('a:\margins.m','1st  loop  VGM  =  [  %4.2f  %4.2f  ]dB\n',vgm2); 
fprintf('a:\margins.m','1st  loop  VPM  =  +  %4.2f  deg\n',vpm2); 
fprintf('a:\margins.m','2nd  loop  VGM  =  [  %4.2f  %4.2f  ]dB\n',vgm3); 
fprintf('a:\margins.m'.'2nd  loop  VPM  =  +  %4.2f  deg\n'.vpm3); 

[agks,bgks,ogks,dgks]=feedbk(agk,-bgk,cgk,-dgk,1); 

[agkt,bgkt,cgkt,dgkt]=feedbk(agk,-bgk,cgk,-dgk,2); 

[svs4,ws4i=sigma(agks,bgks,cgks,dgks); 

[svt4,wt4]=sigma(agkt,bgkt,cgkt,dgkt); 
alphao4=1/max(max(svs4)); 
ao4=1/max(max(svt4)); 
vgms4lo=(1  /(1  +alphao4)); 
vgms4hi=(1/(1-alphao4)); 

if  vgms4hi  <  1e-10;  vgms4hi=1e-10;end; 
vgmt4lo=(1-ao4); 

if  vgmt4lo  <  1e-10;  vgmt4lo=\e-10;end; 
vgmt4hi=(1+ao4); 

vgms4=20*log1 0([vgms4hi  vgms4lo]); 
vgmt4=20*log10([vgmt4hi  vgmt4lo]); 
vpms4=2*asin(alphao4/2)*1 80/pi; 
vpmt4=2'asin(ao4/2)’1 80/pi; 
vpm4=max([vpms4  vpmt4]); 

fprintf('a:\niargins,m','rnv  s-VGM  =  [  %4.2f  %4.2f  ]dB\nrnv  t-VGM  =  [  %4.2f  %4.2f  ]dB\n',vgms4,vgmt4) 
fprintf('a:\margins.m','mv  VPM  =  +  %4.2f  deg\n\n',vpm4) 
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The  following  is  a  MATLAB™  script  file  which  calculates  and  plots  the  singular  values  versus 
frequency  of  GK  and  KG  where  G  is  given  by  [Asys,Busys,Cysys,Dysys]  and  K  is  the  controller 
[Ak,Bk,Ck,Dk]. 


f  igure(fa):subplot(2, 1,1) 

[as,bs,cs,ds]=series(Asys,Busys,Cysys,Dysys,Ak,Bk,Ck,Dk): 

[sv,w]=sigma(as,bs,cs,ds): 

kgsv=20’log10(sv): 

semilogx(w,kgsv,'r'): 

hold:grid;title('SV  plot  of  KG,  H2,  Design  #3r,  RWM‘):xlabel('Frequency  (rad/sec)'),ylabel('Singular  Values  dB') 
axis([5e-3, 1  e3, -50,30]) 

subplot(2,1,2); 

[askg,bskg,cskg,dskg]=serles(Ak,Bk,Ck,Dk,Asys,Busys,Cysys,Dysys); 

[sv,w]=slgma(askg,bskg,cskg,dskg): 

gksv=20*log10(sv); 

semilogx{w,max(gksv),'r'); 

hold;grid;tltle('SV  plot  of  GK,  H2,  Design  #3r,  RWM');xlabel('Frequency  (rad/sec)'),ylabel('Singular  Values  dB') 
axls([5e-3,1e3,-50,30]) 
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theta  -  deg  q  -  deg/sec  alpha  -  deg 


Appendix  C 
Design  Phase  Time  Histories 


Classical  Design 


alpha  vs.  t,  step 


q  vs  t 


■■■'No  Noise 
—  With  Noise 


alpha  vs  t,  IC-5  deg  alpha 
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Figure  C.l. Classical  Design  Time  Histories  (Set  1  of  2) 


C-1 


theta  -  de 


Figure  C.3.  H-2  Design 
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dec  -  deg  gamma  -  deg  Nz  -  g's 
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Nz  vs  t,  Step 
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Figure  C.4.  H-2  Design  1 
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dec  -  deg  gamma  -  deg  Nz  -  g's 


theta  -  deg  q  -  deg/sec  alpha  -  deg 


H-infinity  Design 


"■■No  Noise 


— “  With  Noise 

alpha  vs.  t,  step  alpha  vs  t,  IC-5  deg  alpha 
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Figure  C.5.  H-infinity  Design  Time  Histories  (Set  1  of  2) 
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Figure  C.6.  H-infinity  Design  Time  Histories  (Set  2  of  2) 


Mixed  Design 
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With  Noise 
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Figure  C.7.  Mixed  Design  Time  Histories  (Set  1  of  2) 
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Classical,  H2  ,Hinf,  &  Mixed  -  No  Noise 
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Figure  C.9.  All  Designs  Time  Histories  (Set  1  of  2) 
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Figure  C.IO.  All  Designs  Time  Histories  (Set  2  of  2) 
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Appendix  D 


Data  Acquisition  Parameters 


Parameter 

Units 

Parameter 

Units 

HOURS 

deg 

MINUTES 

deg/s 

SECONDS 

beta_cf 

deg 

MSECONDS 

beta_c 

deg 

hp 

ft 

fes 

lb 

h_dot 

ft/sec 

fas 

lb 

lb 

ftp 

lb 

lb 

des 

in 

fuel_tot 

lb 

das 

deg 

lb 

drp 

in 

oat 

deg 

dec 

deg 

pti 

deg 

dac 

deg 

des_in 

inch 

drc 

deg 

de_virt 

deg 

de 

deg 

dec_fb 

deg 

ds 

deg 

de_test 

deg 

da 

deg 

deg 

dr 

deg 

de_virt2 

deg 

theta 

deg 

dec_m 

deg 

phi 

deg 

decjj 

deg 

psi 

deg 

IBISISIQHIIIIIII 

deg 

P 

q_in 

deg/s 

q 

deg/sec 

kout 

deg 

r 

de.g/sec 

xl 

deg 

nx 

G 

x2 

deg 

ny 

G 

x3 

deg 

nz 

G 

x4 

deg 

nzp 

G 

x5 

deg 

vi 

knots 

x6 

deg 

ve 

knots 

fl6_q 

deg/sec 

vt 

ft/sec 

deg 

lEISS^SIHH 

deg 

1-4 

Table  D.l.  Data  Acquisition  Parameters  Recorded  During  Testing 
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Test  Aircraft:  Lear  25  N1 02VS  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
Data  Reduction:  PC-FRA  Ver  1 .01 
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Data  Analysis  Plots 


Figure  E-1:  Classical  Design  Bode  Plot,  Closed  Loop,  Alpha/PTI  Sweep 


Test  Aircraft:  Lear  25  N102VS  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
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Figure  E-2:  Classical  Design  Bode  Plot,  Closed  Loop,  Pitch  Rate/PTI  Sweep 


Test  Aircraft:  Lear  25  N102VS  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
Data  Reduction:  PC-FRA  Ver  1.01 
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Figure  E-3:  H-2  Design  Bode  Plot,  Closed  Loop,  Alpha/PTI  Sweep 


Test  Aircraft:  Lear  25  N102\/S  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
Data  Reduction:  PC-FRA  Ver  1.01 
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Figure  E-4:  H-2  Design  Bode  Plot,  Closed  Loop,  Pitch  Rate/PTI  Sweep 


Test  Aircraft:  Lear  25  N102VS  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
Data  Reduction:  PC-FRA  Ver  1.01 
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Figure  E-5:  H-Infinity  Design  Bode  Plot,  Closed  Loop,  Alpha/PTI  Sweep 


Test  Aircraft:  Lear  25  N102VS  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
Data  Reduction:  PC-FRA  Ver  1.01 


Figure  E-^ 


Test  Aircraft:  Lear  25  N102VS  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input;  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
Data  Reduction:  PC-FRA  Ver  1.01 
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Figure  E-7:  Mixed  Design  Bode  Plot,  Closed  Loop,  Alpha/PTI  Sweep 


Test  Aircraft:  Lear  25  N102\/S  Dates:  30  Sept  -  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  PTI  Frequency  Sweep 

Configuration:  Gear-Down;  Flaps-20%  0.1-10  Hz;  Amplitude  0.25  Degrees 
Data  Reduction:  PC-FRA  Ver  1.01 
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Figure  E-8:  Mixed  Design  Bode  Plot,  Closed  Loop,  Pitch  Rate/PTI  Sweep 


Figure  E-9:  Step  Input  Time  Histories,  Classical  Design 


E-10 


Figure  E-11:  Step  Input  Time  Histories,  H-Infinity  Design 


Test  Aircraft:  Lear  25  N102VS  Date:  3  Oct  96 

Altitude:  4,000  -  6,000  ft  MSL  Input:  Step,  1"  stick  displacement 

Configuration:  Gear  -  Down;  Flaps  -  20  %  - Flight  Test  Data 


MATLAB  Model  Prediction 


1  2  3  Time  (sec)  4  5  6 


Figure  E-12:  Step  Input  Time  Histories,  Mixed  Design 
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APPENDIX  F 


Test  Card  and  Rating  Scales 


HAVE  INFINITY  TEST  CARD 


^AL  PILOT: 


TURB: 


APP  SPEED:  I  D  SPEED: 

_ /  / 


TASK:  Intercept  ILS  glideslope  and  offset  300’  laterally.  At  200’  AGL,  correct 
to  centerline  by  100’  AGL  using  30  to  45  degrees  of  bank.  Intercept  a  visual 
elidepath  to  touchdown  in  desired  box. 


CH:  PIO: 


LANDING: 

Desired 

INITIAL  RESPONSE: 

SLOW 

1 _ 

SENSITIVITY: 

LOW 

I _ 

PREDICTABILITY: 

BAD 

1 _ 

LOW 

WORKLOAD:  I _ 

COMPENSATION  (Type): 

UNDESIRABLE  MOTIONS  (WhenH’ype): 

PROBLEMS  EFFECTING  TASK: 

Wind 

Turbulence 

Offset 

Glideslope 

Correction 

Lateral  Control 

Flare 

Power  Control 
Touchdown 
VSS  Trip 

GEN  HANDLING  QUALITIES: 


Figure  FI  Sample  Test  Card 
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COOPER-HARPER  RATING  SCALE 


AIRCRAFF  DEMANDS  ON  THE  PILOT  IN  SELEa’ED  PlI 
.RACTERISTICS  TaSKOR  REQi;iRED  OPEBATON  BA' 
Excellent  Pilot  compensation  not  a  tactor  for  B 

Hiithiy  desiaMn’f  ilcsircJ  penorraince  m 

Good  Pilot  compensation  not  a  tactor  tor  B 

Neglibable  deficiencies  desired  performance  B 

Fair-  Some  inildK  Minimal  pilot  Lompensaiion  require  J  B 

unpleasant  deficiencies  for  desired  performance  B 


Miimr  but  .innovine  Desired  |)enom'-ince  requires 

deficiencies  moderate  pilot  compensation 

Moderately  objectionable  Adequate  performance  requires 

deficencies  considerable  pilot  compensation 

Very  objectionable  but  Adequate  performance  requires 

tolerable  deficiencies  extensive  pilot  con^nsation 


Deficiencies 

warrant 

inprovement 


Adequate  performance  not  attainable  wit 
maximum  tolerable  pilot  compensation. 
Controllability  not  in  question 
Considerable  pilot  compensation  is 
required  for  control 
Intense  pilot  compensation  is  required 
to  retain  control 


adcqujieV. 
lerforamnce  attaina^ 
with  a  tolerable  pilo^ 
workload 


Deficiencies 

require 

improvement 


Major 

deficencies 

Major 

deficencies 


ViQ  ^  Improvement  g' §Major 

mandatory  |!^deficencies 


Is  It 

Controllable? 


Control  will  be  lost  during  some 
portion  of  required  operation 


Pilot  Decisons 


Figure  F2  Cooper  Harper  Rating  Scale 
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No 


2 


PIO  1  -  No  tendency  for  pilot  to  induce  undesirable  motion. 

PIO  2  -  Undesirable  motion  tends  to  occur  when  pilot  initiates  abrupt  maneuvers  or  attempts  tight  control. 
These  motions  can  be  prevented  or  eliminated  by  pilot  technique. 

PIO  3  -  Undesirable  motions  easily  induced  when  pilot  initiates  abrupt  maneuvers  or  attempts  tight  control. 
These  motions  can  be  prevented  or  eliminated,  but  only  at  sacrifice  to  task  performance  or  through 
considerable  pilot  attention  and  effort. 

PIO  4  -  Oscillations  tend  to  develop  when  pilot  initiates  abrupt  maneuvers  or  attempts  tight  control.  Pilot 
must  reduce  gain  or  abandon  task  to  recover. 

PIO  5  -  Divergent  oscillations  tend  to  develop  when  pilot  initiates  abrupt  maneuvers  or  attempts  tight 
control.  Pilot  must  reduce  gain  or  abandon  task  to  recover. 

PIO  6  -  Disturbance  or  normal  control  may  cause  divergent  oscillation.  Pilot  must  open  control  loop  by 
releasing  or  freezing  the  stick. 
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Appendix  G 


Pilot  Comments  Data  Base 


CONFIGURATION  ID  Classical 


Mission  date:  30  Sep  96  Eval  pilot:  Fittante  Flight  #  1 

Handling  qualities:  The  aircraft  was  responsive  and  predictable.  Overall  the  pilot  stated  that  this 
configuration  had  nice  handling  qualities. 

Landing:  The  nose  tracked  well  and  the  task  was  easy  to  perform. 


Appr 

Landing  zone 

Turb 

Wind 

C-H  Rating 

PIO  Rating 

1 

Desired 

none 

calm 

3 

1 

2 

Desired 

none 

calm 

3 

1 

3 

Desired 

none 

calm 

3 

1 

Notes  on  C-H:  All  four  had  desired  performance  and  minimal  workload. 
Notes  on  PIO  Rating:  No  undesirable  motions  were  noted. 


Mission  date:  1  Oct  96  Eval  pilot:  Edwards  Flight  #  2 

Handling  qualities:  The  aircraft  had  a  quick  response  and  was  fairly  predictable  with  moderate  to  high 
sensitivity. 

Landing:  There  was  no  tendency  to  balloon  or  bounce  with  this  configuration. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  Light  vrb/4  4  2 

2  Desired  Light  vrb/4  4  1 

Notes  on  C-H:  Moderate  workload  was  required  because  the  response  of  this  configuration  required  the 
pilot  to  make  quick,  abrupt  inputs  to  maintain  glide  path. 

Notes  on  PIO  Rating:  There  was  a  tendency  to  over  flare  for  the  first  approach  and  no  unwanted  motions 
were  perceived  on  the  second  approach. 
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CONFIGURATION  ID  Classical 


Mission  date:  2  Oct  96  Eval  pilot:  Malacrida  Flight  #  3 

Handling  qualities:  The  aircraft  predictability  was  good  and  the  sensitivity  was  nice.  The  stick  was  a  bit 
heavy. 

Landing:  No  problems  noted  on  landing. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  180/3  4  2 

2  Desired  none  180/3  4  2 

Notes  on  C-H:  The  pilot  described  his  workload  as  moderate,  but  also  commented  that  this  is  a  borderline 
level  one  design. 

Notes  on  PIO  Rating:  A  slight  bobble  resulted  during  the  ILS  glideslope  capture  and  resulted  in  a  PIO  2. 

Mission  date:  2  Oct  96  Eval  pilot:  Fittante  Flight  #  4 

Handling  qualities:  The  aircraft  had  good  response  with  fairly  light  stick  forces. 

Landing:  The  aircraft  nose  stayed  down  after  landing. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  vrb/5  3  1 

Notes  on  C-H:  Minimal  workload  was  required. 

Notes  on  PIO  Rating:  No  unwanted  motions  were  observed. 


Mission  date:  3  Oct  96  Eval  pilot:  Bouchard  Flight  #  5 

Handling  qualities:  The  aircraft  had  moderate  to  slow  initial  response,  with  moderate  sensitivity, 
predictability  and  workload. 

Landing:  The  pilot  felt  he  flared  a  little  to  high  which  resulted  in  a  firm  landing  for  the  first  approach  and 
the  aircraft  had  excessive  aft  stick  required  for  flare. 

Appr  Landing  zone  Turb  Wind  C-H  Rating 

1  Adequate  none  220/6  5 

2  Desired  none  220/6  4 

Notes  on  C-H:  A  CH  5  was  given  for  the  first  approach  because  of  adequate  performance,  but  the  pilot 
commented  that  the  handling  qualities  were  more  representative  of  CH  3.  The  second  approach  was  given  a 
CH  4  because  of  moderate  compensation  in  the  form  of  excessive  aft  stick  required  to  flare. 
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CONFIGURATION  ID  Hz 


Mission  date:  30  Sep  96  Eval  pilot:  Fittante  Flight  #  1 

Handling  qualities:  The  aircraft  had  heavy  stick  forces  and  didn’t  want  to  move.  It  was  very  stable  with  a 
dead  beat  response,  and  was  unresponsive.  Sharp  inputs  caused  a  slight  pitch  bobble. 

Landing:  The  aircraft  was  not  evaluated  for  control  after  touchdown  since  the 
evaluation  pilot  tripped  off  the  VSS. 


Appr 

Landing  zone 

Turb 

Wind 

C-H  Rating 

PIO  Rating 

1 

Desired 

none 

calm 

4 

2 

2 

Desired 

none 

calm 

4 

2 

3 

Desired 

none 

calm 

4 

2 

Notes  on  C-H:  On  the  second  approach  the  initial 

correction 

was  input  with  the  aircraft  one  dot  above 

glideslope  causing  a  third  approach. 

Notes  on  PIO  Rating:  The  rating  of  2  was  given  due  to  the  slight  pitch  bobble. 


Mission  date:  1  Oct  96  Eval  pilot:  Edwards  Flight  #  2 

Handling  qualities:  The  aircraft  was  very  quick  in  response,  jittery,  and  required  low  gain  smooth  inputs 
to  obtain  the  desired  results. 

Landing:  During  the  flare  the  stick  had  to  be  rapidly  brought  forward  to  lower  the  nose. 


Appr 

Landing  zone 

Turb 

Wind 

C-H  Rating 

PIO  Rating 

1 

Desired 

Light 

160/3 

5 

3 

2 

Desired 

Light 

230/6 

4 

1 

Notes  on  C-H:  Sensitivity  to  control  inputs  made  the  aircraft  unpredictable  which  increased  workload  for 
the  first  approach.  For  the  second  approach  workload  decreased  due  to  a  learning  curve. 

Notes  on  PIO  Rating:  Undesirable  motions  caused  the  rating  of  3  for  the  first  approach,  but  were  not 
perceived  on  the  second  approach  which  was  given  a  rating  of  1, 
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CONFIGURATION  ID  Hj 


Mission  date:  2  Oct  96  Eval  pilot:  Malacrida  Flight  #  3 

Handling  qualities:  The  aircraft  had  poor  pitch  predictability  for  higher  gain  inputs.  The  initial  response 
of  the  aircraft  was  slow  and  sensitivity  was  low. 

Landing:  The  pilot  needed  small  pulse  inputs  to  flare.  The  aircraft  was  barely  controllable  on  the  ground 
requiring  high  pilot  compensation. 


Appr 

Landing  zone 

Turb 

Wind 

C-H  Rating 

PIO  Rating 

1 

Desired 

none 

280/6 

8 

4 

2 

Desired 

none 

280/6 

8 

4 

Notes  on  C-H:  The  pilot  had  full  forward  stick  applied  at  touchdown  and  still  had  difficulty  keeping  the 
aircraft  on  the  ground. 

Notes  on  PIO  Rating:  Significant  non  oscillatory,  undesirable  motions  were  seen  on  the  ground. 


Mission  date:  2  Oct  96  Eval  pilot:  Fittante  Flight  #  4 


Handling  qualities:  The  aircraft  had  heavy  stick  forces  and  was  sluggish.  The  heavier  forces  and  delayed 
response  produced  a  little  oscillation  on  final. 


Landing:  At  approximately  two  seconds  after  touchdown  full  forward  stick  was 

required  to  keep  the  nose 

from  pitching  up. 

Appr 

Landing  zone 

Turb 

Wind  C-H  Rating 

PIO  Rating 

1 

Desired 

none 

vrb/5  6 

3 

2 

Desired 

none 

vrb/5  8 

3 

Notes  on  C-H:  Compensation  was  mainly  required  for  the  sluggish  pitch  response  and  the  pitch  up  after 
landing.  On  the  second  approach  controllability  was  in  question  after  touchdown. 

Notes  on  PIO  Rating:  The  rating  of  3  was  given  for  the  undesired  motion  on  final. 


Mission  date:  3  Oct  96  Eval  pilot:  Bouchard  Flight  #  5 

Handling  qualities:  The  aircraft  had  low  sensitivity,  bad  predictability  and  high  workload.  On  final 
approach  the  aircraft  seemed  to  handle  fine. 

Landing:  On  the  straight  in  at  touchdown  the  aircraft  bounced  and  began  an  uncontrollable  pitch  up 
maneuver  upon  which  the  safety  pilot  disengaged  the  VSS  and  initiated  a  go  around. 

Notes  on  C-H:  Although  this  test  did  not  consider  the  straight  in  an  operational  task  with  regards  to 
Cooper  Harper  ratings,  the  configuration  was  deemed  uncontrollable. 
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CONFIGURATION  ID  Hz 


Mission  date:  3  Oct  96  Eval  pilot:  Malacrida  Flight  #  6 

Handling  qualities:  The  initial  response  of  the  aircraft  was  good  but  the  sensitivity  was  low  to  moderate 
and  the  predictability  was  bad. 

Landing:  At  approximately  two  seconds  after  touchdown  full  forward  stick  was  required  to  keep  the  nose 
from  pitching  up. 


Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Not  adequate  none  vrb/5  7  4 

Notes  on  C-H:  The  CH  7  was  based  on  the  long  landing  and  the  problems  in  the  flare. 

Notes  on  PIO  Rating:  The  task  was  abandon  due  to  pitch  bobble  and  uncommanded  pitch  up. 


Mission  date:  3  Oct  96  Eval  pilot:  Bouchard  Flight  #  6 

Handling  qualities:  The  aircraft  had  good  initial  response  with  slightly  low  sensitivity  and  bad 
predictability. 

Landing:  Forward  stick  required  in  flare  and  landing  to  keep  the  aircraft  from  ballooning. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  vrb/5  10  1 

2  Desired  none  vrb/5  9  1 

Notes  on  C-H:  The  aircraft  was  uncontrollable  after  touch  down  for  the  first  approach  and  with  full  stick 
on  the  second  approach  the  aircraft  marginal  control  on  roll  out. 
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CONFIGURATION  ID  E, 


Mission  date:  1  Oct  96  Eval  pilot:  Edwards  Flight  #  2 

Handling  qualities:  The  stick  inputs  were  sensitive  with  a  quick  initial  response,  glidepath  was  not 
difficult  to  maintain 

Landing:  The  aircraft  wanted  to  float  in  the  flare.  During  the  flare  the  stick  had  to  be  brought  forward  to 
prevent  ballooning. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Safety  trip  Light  190/5  7  2 

2  Adequate  Light  190/5  7  2 

Notes  on  C-H:  On  the  first  rated  approach  the  safety  pilot  took  control  prior  to  touchdown  due  to  a  high 
sink  rate.  On  the  second  approach  touched  down  in  the  adequate  box  with  the  aircraft  ballooning. 

Notes  on  PIO  Rating:  A  PIO  of  2  was  given  due  to  unwanted  motion  in  the  flare. 


Mission  date:  2  Oct  96  Eval  pilot:  Malacrida  Flight  #  3 

Handling  qualities:  The  aircraft  was  predictable  and  had  good  pitch  response  until  the  pilot  gains  were 
increased. 

Landing:  The  aircraft  had  a  tendency  to  float  on  the  runway.  Too  much  forward  stick  was  required  during 
the  flare. 


Appr 

Landing  zone 

Turb 

Wind 

C-H  Rating 

PIO  Rating 

1 

Not  Adequate 

none 

280/6 

1 

3 

2 

Adequate 

none 

280/6 

5 

3 

3 

Not  Adequate 

none 

280/6 

7 

3 

Notes  on  C-H:  Workload  was  not  an  issue,  poor  performance  drove  the  CH  5  and  7’s. 

Notes  on  PIO  Rating:  The  unwanted  motion  of  the  heave  and  float  at  flare  resulted  in  a  PIO  3. 
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CONFIGURATION  ID  H. 


Mission  date:  2  Oct  96  Eval  pilot:  Fittante  Flight  #4 

Handling  qualities:  The  aircraft  had  good  aircraft  response  through  the  correction  to  centerline  with  light 
pitch  forces. 

Landing:  The  aircraft  wanted  to  float  and  forward  stick  was  required  for  touchdown.  Full  forward  stick 
was  required  to  keep  the  nose  on  the  runway. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  vrb/5  8  3 

2  Desired  none  vrb/5  8  3 

Notes  on  C-H:  The  aircraft  was  just  barely  controllable  for  both  landings  after  touchdown. 

Notes  on  PIO  Rating:  A  PIO  rating  of  3  was  assigned  for  the  easily  induced  jittery  aircraft  response 
during  the  flare  as  stick  forces  became  lighter. 


Mission  date:  2  Oct  96  Eval  pilot:  Edwards  Flight  #  4 

Handling  qualities:  The  lateral  correction  was  easily  accomplished  with  heavy  stick  forces. 

Landing:  The  aircraft  had  a  tendency  to  balloon  in  the  flare  and  full  forward  stick  was  required  to  keep  the 
aircraft  on  the  ground. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  200/5  8  2 

2  Desired  none  200/5  7  2 

Notes  on  C-H:  Pilot  workload  to  prevent  the  aircraft  from  pitching  up  drove  the  CH  ratings. 

Notes  on  PIO  Rating:  A  PIO  of  2  was  given  due  to  unwanted  motion  during  and  after  touchdown. 


Mission  date:  3  Oct  96  Eval  pilot:  Bouchard  Flight  #  5 

Handling  qualities:  The  aircraft  had  moderate  to  quick  initial  response,  with  moderate  to  highly  sensitive, 
slightly  unpredictable  and  moderate  workload. 

Landing:  Forward  stick  was  required  to  keep  the  nose  down. 

Appr  Landing  zone  Turb  Wind  C-H  Rating 

1  Desired  none  220/6  7 

2  Desired  none  220/6  9 

Notes  on  C-H:  The  design  was  considered  level  3  because  of  the  roll  out  characteristics.  The  first  received 
a  CH  7  because  of  extensive  pilot  compensation  without  questionable  controllability.  The  second  received 
a  CH  9  because  controllability  was  in  question  after  touchdown. 


G-7 


CONFIGURATION  ID  EL 


Mission  date:  3  Oct  96  Eval  pilot:  Malacrida  Flight  #  6 

Handling  qualities:  The  aircraft  had  adequate  initial  response  and  sensitivity  but  it  had  bad  predictability. 

Landing:  A  forward  stick  was  required  at  flare  to  keep  the  aircraft  from  ballooning. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Not  adequate  none  vrb/5  7  3 

Notes  on  C-H:  The  aircraft  landed  long. 

Notes  on  PIO  Rating:  The  radng  of  3  was  given  for  the  undesired  motion  on  final. 

Mission  date:  3  Oct  96  Eval  pilot:  Bouchard  Flight  #  6 

Handling  qualities:  The  aircraft  had  adequate  initial  response  with  low  sensitivity  and  bad  predictability. 

Landing:  After  touchdown  forward  stick  was  required  to  keep  the  nose  from  pitching  up. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Adequate  none  vrb/5  6  1 

Notes  on  C-H:  The  aircraft  landed  in  the  adequate  zone  with  high  workload. 
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CONFIGURATION  ID  Mixed  Hz/BL 


Mission  date:  1  Oct  96  Eval  pilot:  Edwards  Flight  #  2 

Handling  qualities:  The  aircraft  handled  well  coming  down  final  and  during  the  offset  task.  The  aircraft 
initial  response  seemed  slow. 

Landing:  During  the  flare  the  stick  had  to  be  brought  forward  to  prevent  ballooning. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  230/6  5  2 

2  Desired  none  230/6  8  2 

Notes  on  C-H:  On  the  first  approach  a  CH  5  was  given  due  to  high  workload.  On  the  second  approach  a 
CH  8  was  given  due  to  unpredictable  aircraft  motion  after  touchdown 

Notes  on  PIO  Rating:  A  PIO  of  2  was  given  due  to  unwanted  motion  after  touchdown. 


Mission  date:  2  Oct  96  Eval  pilot:  Malacrida  Flight  #  3 

Handling  qualities:  The  aircraft  felt  heavy  with  a  sluggish  response  and  was  not  terribly  predictable. 

Landing:  The  aircraft  exhibited  a  heave  with  the  flare. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Adequate  none  280/6  7  3 

2  Not  adequate  none  280/6  7  3 

Notes  on  C-H:  The  pitot  gave  CH  7  because  of  the  high  workload  and  extensive  compensation.  The  pilot 
had  to  consciously  lower  his  gains  and  apply  his  inputs  early. 

Notes  on  PIO  Rating:  The  unwanted  motion  of  the  heave  at  flare  resulted  in  a  PIO  3. 


Mission  date:  2  Oct  96  Eval  pilot:  Fittante  Flight  #  4 

Handling  qualities:  The  aircraft  had  slightly  high  stick  forces  and  pitch  control  was  solid  during  transition 
to  landing  until  just  prior  to  touchdown. 

Landing:  At  touchdown  there  was  a  sudden  reversal  in  pitch  force,  requiring  full  forward  stick  to  keep  the 
nose  on  the  runway. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  vrb/5  8  1 

2  Desired  none  vrb/5  8  1 

Notes  on  C-H:  CH  8  was  given  for  controllability  during  and  after  touchdown. 

Notes  on  PIO  Rating:  The  PIO  rating  of  1  was  given  because  there  was  no  longer  a  jittery  response  during 
the  transition  to  landing. 
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CONFIGURATION  ID  Mixed  Hi/H^ 


Mission  date:  2  Oct  96  Eval  pilot:  Edwards  Flight  #  4 

Handling  qualities:  The  initial  response  was  good  with  adequate  sensitivity  and  bad  predictability. 

Landing:  Full  forward  stick  was  applied  after  touchdown,  but  the  aircraft  became  airborne  again  after  both 
landings. 


Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Desired  none  vrb/5  8  2 

2  Desired  none  vrb/5  8  2 

Notes  on  C-H:  Considerable  pilot  compensation  was  required. 

Notes  on  PIO  Rating:  A  PIO  of  2  was  given  due  to  unwanted  motion  after  touchdown. 


Mission  date:  3  Oct  96  Eval  pilot:  Bouchard  Flight  #  5 

Handling  qualities:  The  aircraft  had  moderate  initial  response,  with  low  to  moderate  sensitivity,  bad 
predictability  and  high  workload. 

Landing:  At  touchdown,  during  the  straight  in,  the  aircraft  bounced  and  began  an  uncontrollable  pitch  up 
maneuver  similar  to  configuration  L.  At  this  time  the  pilot  requested  the  safety  pilot  to  take  control  of  the 
aircraft. 

Notes  on  C-H:  Although  this  test  did  not  consider  the  straight  in  an  operational  task  with  regards  to  Cooper 
Harper  ratings,  the  configuration  was  deemed  uncontrollable. 


Mission  date:  3  Oct  96  Eval  pilot:  Malacrida  Flight  #  6 

Handling  qualities:  The  aircraft  had  adequate  initial  response  but  moderately  low  sensitivity  and  bad 
predictability. 

Landing:  After  touchdown  forward  stick  was  required  to  keep  the  nose  from  pitching  up. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Not  adequate  none  vrb/5  7  3 

Notes  on  C-H:  The  aircraft  landed  long. 

Notes  on  PIO  Rating:  The  rating  of  3  was  given  for  the  undesired  motion  on  final. 
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CONFIGURATION  ID  Mixed  Hj/H^ 


Mission  date:  3  Oct  96  Eval  pilot:  Bouchard  Flight  #  6 

Handling  qualities:  The  aircraft  had  adequate  initial  response  and  sensitivity  and  bad  predictability. 

Landing:  After  touchdown  full  forward  stick  was  required  to  keep  the  nose  from  pitching  up. 

Appr  Landing  zone  Turb  Wind  C-H  Rating  PIO  Rating 

1  Adequate  none  vrb/5  6  1 

Notes  on  C-H:  Landed  in  the  adequate  zone  and  had  high  workload. 
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Number  of  Ratings 


Appendix  H 


PIO  Rating 


Figure  H.2.  Classical  Flight  Control  Design  Pilot  Induced  Oscillation  Ratings 
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Number  of  Ratings  Number  of  Ratings 


H2  Design: 


Cooper  Harper  Ratings 

Figure  H.3.  H2  Flight  Control  Design  Handling  Quality  Ratings 


PIO  Rating 

Figure  H.4.  H2  Flight  Control  Design  Pilot  Induced  Oscillation  Ratings 
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Number  of  Ratings  Number  of  Ratings 


H>,  Design: 


Test  Aircraft:  Lear  25,  N102VS  Dates:  30  Sept  -  3  Oct  1996 

Configuration:  Gear  Down,  Raps  20%  Data  Source:  Rioht  Test 


1  23456789  10 


Cooper  Harper  Ratings 

Figure  H.5.  H<x.  Flight  Control  Design  Handling  Quality  Ratings 


PIO  Rating 


Figure  H.6.  Hco  Flight  Control  Design  Pilot  Induced  Oscillation  Ratings 


Number  of  Ratings  Number  of  Ratings 


Mixed  Design: 


Cooper  Harper  Rating 

Figure  H.7.  Mixed  Flight  Control  Design  Handling  Quality  Ratings 


Figure  H.8.  Mixed  Flight  Control  Design  Pilot  Induced  Oscillation  Ratings 
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